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ABSTRACT 


The  PHALANX  Close-In  Weapons  System  (CIWS)  is  a  very  important  part  of  the  anti¬ 
ship  missile  defense  structure  of  the  U.S.  Navy.  The  PHALANX  gun  system  currently 
experiences  random  and  variable  bullet  dispersion  which  diminishes  its  ability  to  destroy  its 
intended  targets.  The  more  that  is  understood  about  parameters  that  control  PHALANX 
gun  dispersion  performance,  the  more  the  Navy  can  increase  its  performance  and  hence,  ship 
survivability. 

The  research  described  in  this  thesis  is  an  attempt  to  correlate  measured  bullet 
dispersion  with  transverse  barrel  tip  displacement  data  derived  from  measured  barrel  tip 
acceleration  data  taken  on  a  firing  PHALANX  gun  system.  The  goal  of  this  project  is  to 
relate  barrel  tip  vibration,  which  can  be  measured  on  a  gun  at  any  location,  to  bullet 
dispersion,  which  can  only  be  measured  on  a  test  range. 

Acceleration  measurements  were  made  on  a  firing  PHALANX  gun  at  the  Navy  Air 
Weapons  Station,  China  Lake,  CA.  From  these  measurements  were  derived  the  vertical  and 
horizontal  components  of  gun  barrel  tip  displacement  by  direct  time  integration  and  by  a 
fourier  transform  procedure.  It  was  found  that  the  resulting  displacement  time  records  were 
contaminated  by  a  large,  spurious,  low  frequency  contribution,  the  source  of  which  is 
unknown.  It  is  speculated  that  it  is  most  likely  due  to  the  effect  of  temperature  on  the 
accelerometers.  A  procedure  was  developed  to  remove  the  low  frequency  trend  from  the 
displacement  time  records  using  cubic  spline  smoothing.  The  resulting  horizontal  and 
vertical  barrel  tip  displacement  records  were  compared  to  bullet  dispersion  data  gathered 
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during  the  gun  firing.  Although  there  was  a  good  correlation  between  the  time  of  firing 
obtained  from  the  dispersion  data  and  the  time  of  maximum  vertical  displacement, 
correlation  between  dispersion  data  and  extracted  lateral  barrel  tip  displacement  was 
disappointing.  Suggestions  are  made  to  improve  the  procedure  for  extracting  the  barrel  tip 
lateral  displacement  from  vibration  measurements  on  a  firing  gun. 
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I.  INTRODUCTION 


All  gun/ammunition  systems  have  a  dispersion  figure.  Dispersion  refers  to  the 
lateral  rms  error  in  putting  rounds  on  a  target,  usually  measured  as  an  angle.  If  a  number 
of  rounds  are  fired  from  a  weapon  against  a  target  there  will  be  dispersion  in  both  the 
horizontal  (X)  and  vertical  (Y)  lateral  directions.  For  small  caliber  weapons,  lateral 
dispersion  tends  to  be  small  and  equal  in  both  directions.’  It  is  known  that  the  horizontal 
and  vertical  dispersions  are  normally  distributed  if  the  X  and  Y  dispersion  figures  are 
equal.  ^  Typically,  dispersion  is  measured  in  terms  of  the  standard  deviation  of  a 
Gaussian  distribution.  It  can  also  be  expressed  in  terms  of  the  radius  of  a  circle  which 
encompasses  80  percent  of  the  rounds  in  a  pattern. 

There  is  a  considerable  difference  in  hit-probability  between  firing  a  new  weapon 
and  one  with  wear,  between  firing  amour-piercing  and  high-explosive  ammunition,  and 
between  firing  single-shot  and  in  bursts.  When  a  weapon  is  fired  in  bursts,  as  opposed  to 
single  shots,  dispersion  increases  by  50  per  cent.^  For  the  PHALANX  gun  system  the 
dispersion  is  quoted  to  be  1.4  mradians  at  one  sigma.  [Ref  1] 

Dispersion  data  for  a  firing  PHALANX  gun  system  were  collected  at  the  Naval 
Air  Weapon  Station,  China  Lake,  CA.  The  data  were  collected  in  conjunction  with 
accelerometer  data.  The  intent  of  this  research  is  to  correlate  bullet  dispersion  and 
transverse  barrel  tip  displacement  data  obtained  on  a  firing  PHALANX  gun 
system. 


’  International  Defense  Review  1 1/91 
^  International  Defense  Review  11/91 
^  International  Defense  Review  11/91 
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A.  PHALANX  CLOSE-IN  WEAPON  SYSTEM 

The  PHALANX  Close-In  Weapon  System  (CIWS)  program  was  initiated  in  the 
late  1960's  to  provide  the  Navy  with  close-in  autonomous  defense  against  low-altitude 
Anti-Ship  Missiles  (ASMs).  Low-altitude  ASMs  were  considered  the  primary  threat  of 
the  1970s  because  they  could  penetrate  the  other  principal  defense  systems  in  the  fleet. 
Since  then,  there  have  been  significant  changes  to  the  ASM  threat. 

PHALANX  was  first  introduced  to  the  fleet  as  a  Block  0  system.  Since  then. 
Block  0  systems  have  formed  the  backbone  of  the  Navy’s  current  CIWS  defense 
program.  As  the  program  progresses,  so  does  CIWS  technology.  The  PHALANX 
program  has  introduced  system  improvements  through  upgrades,  modifications,  and  new 
subsystems  supported  by  an  evolving  test  and  integration  schedule.  Improvements  from 
Block  0  to  Block  1  Baseline  2  include  pneumatic  gun  drive  to  allow  higher  firing  rates 
with  higher  reliability,  improved  search  sensitivity  through  new  standing  wave  planar 
array  search  antennas,  and  the  gun  muzzle  restraint  to  control  gun  dispersion.  [Ref  2] 

1.  System  Description 

The  PHALANX  gun  system  is  designed  as  a  stand-alone  integrated  weapon 
system  that  encompasses  functions  usually  performed  by  separate,  independent  systems. 
As  an  integrated  package,  it  carries  out  search,  detection,  target  declaration,  tracking, 
threat  elevation,  declaration  of  engagement,  fire  control,  and  target  kill  assessment.  The 
integrated  concept  (shown  in  Figure  1-1)  provides  the  fast  reaction  time  required  for  the 
CIWS  mission. 

PHALANX  consists  of  six  major  assemblies  and  components.  Five  of  these  are 
the  radar  and  servo  assembly,  gun  assembly,  mount  and  train  drive  platform,  barbette 
equipment  assembly,  and  electronics  enclosure.  The  sixth  component  is  made  up  of  the 
local  control  panel  and  the  remote  control  panel. 

a.  On-Mount  Fire  Control  Assembly 

The  on-mount  fire  control  assembly  performs  the  following  functions: 
target  search,  detection  and  threat  declaration,  track  acquisition,  target  track  and 
measurement  of  range,  velocity  and  angle,  target  prediction,  lead  angle  computation,  gun 
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aiming  and  firing,  projectile  detection,  measurement  of  projectile  velocity  and  angle,  and 
gun  aim  bias  correction. 


AMCMILV 


Figure  1-1  PHALANX  System  Components 

Fast  reaction  times  are  achieved  by  using  search  and  track  radars  along  with  the 
weapon  control  computer  to  achieve  target  acquisition.  This  is  accomplished  by 
maintaining  the  search  platform  in  a  horizontally-stabilized  position  and  salved  to  a 
vertical  gyro  so  that  targets  may  be  sorted  and  correlated  according  to  range,  range  rate, 
and  angular  position.  The  track  antenna  reference  system,  with  its  own  rate-integrating 
gyro  package,  is  space-stabilized  upon  command  of  the  weapon  control  computer  to 
ensure  fast  handover  and  track  acquisition. 

The  search  antenna  subsystem  consists  of  standing  wave  antennas  mounted  to  the 
search  platform  at  specific  elevation  angles  to  provide  a  large  elevation  coverage.  The 
combination  of  the  track  radar,  antenna  control,  fire  control  computer  solution  and  mount 
servo  defines  the  engagement  control  system  for  the  PHALANX  system.  Using  the 
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antenna  position  as  a  reference,  the  difference  between  gun  position  and  director  solution 
provides  an  error  signal  that  generates  a  gun  pointing  correction.  This  technique 
eliminates  the  need  for  independent  antenna  position  measurement  and  eliminates  deck 
reference  with  its  flexure  problems.  This  technique  allows  extremely  accurate  angle  error 
measurement  between  mount  and  antenna. 

b.  Gun  Assembly 

The  heart  of  the  PHALANX  weapon  system  is  the  versatile  M61 A1 
gatling  gun,  providing  a  high  rate  of  fire  (variable  3000  or  4500  rounds  per  minute)  with 
a  specifically-designed  high  kinetic  energy  round,  to  provide  the  necessary  lethality.  The 
six-barrel  Gatling  gun  was  adopted  from  the  Air  Force  M61  Vulcan  gun  series  used  in 
several  types  of  aircraft  and  ground  mounted  for  airfield  defense.  '*  The  gun  assembly 
(Figure  1-2)  is  made  up  of  five  major  sub-assemblies.  The  gun  is  electronically 
controlled  and  pneumatically  driven.  The  gun  consists  of  a  rotating  cluster  of  six  barrels 
with  a  breech  bolt  for  each  barrel.  The  ammunition  subsystem  includes  all  components 
that  store,  feed  and  transport  rounds  and  empty  cartridge  cases  within  the  system. 

c.  Mount  and  Train-Drive  Assembly 

The  mount  provides  structural  support  to  the  gun  assembly  and  the  radar 
servo  assembly  and  houses  the  servo  drive  components  that  provide  elevation  control  of 
the  gun.  The  train  drive  assembly  is  a  cast  structure  that  contains  the  train  bearing.  This 
assembly  is  a  combination  gear  and  bearing  which  serves  as  the  interface  between  the  gun 
mount  and  barbette.  The  train  drive  gear  box  is  mounted  to  the  platform  assembly  and 
provides  azimuth  drive  to  the  gun  mount  through  the  train  bearing.  An  instrument 
gearbox  is  also  attached  to  the  train  drive  platform  and  provides  position  feedback  to  the 
azimuth  servo  loop  system  and  fire  interrupt  for  no  fire  sectors. 


*  The  Ships  and  Aircraft  of  the  U.S.  Fleet,  Thirteenth  edition. 
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Figure  1-2  Gun  Assembly 


d.  Barbette 

The  barbette  houses  the  transmitter,  transmitter  power  supply,  pneumatic 
drive  power  supply,  pneumatic  electronic  control  unit,  power  controller  transformer 
assembly,  and  the  electronic  control  group  (ECG).  The  PHALANX  CIWS  provides 
environmental  control  for  its  weapon  group  equipment,  receiving  only  seawater  and 
electrical  power  from  the  ship  for  this  function.  The  ECG  provides  air  to  the  system 
electronics  for  temperature  and  condensation  control,  dry  air  for  waveguide  pressurization 
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and  a  water/glycol  mixture  for  transmitter  cooling. 

e.  Electronics  Enclosure 

The  electronics  enclosure  provides  weather  and  electromagnetic 
interference  protection  for  the  rest  of  the  electronics  that  includes  the  signal  processor, 
signal  generator,  gun  control  unit,  mount  elevation  and  train  power  drives,  search  and 
track  servo  electronics,  power  supply  control  group  and  the  weapon  control  group 
computer. 

/.  CIWS  MK 149  Round 

The  PHALANX  CIWS  utilizes  20  millimeter  projectiles  that  prove  to  be 
lethal  against  an  incoming  missile.  This  projectile  must  have  a  high  muzzle  velocity  and 
be  capable  of  penetrating  numerous  components  in  a  missile  forebody.  Having 
penetrated  or  otherwise  exposed  the  warhead,  it  must  impact  the  explosive  material  with 
the  correct  relationship  of  residual  velocity  and  cross-sectional  area  but  also  cause  a 
warhead  reaction.  Development  tests  have  resulted  in  the  MK149  armor  piercing 
discarding  sabot  round  (Figure  1-3),  which  is  a  sub-caliber,  spin-stabilized  tungsten 
penetrator.  At  firing,  the  pusher  imparts  spin  to  the  projectile  and  forms  a  gas  seal.  After 
leaving  the  muzzle,  both  the  lightweight  pusher  and  sabot  are  discarded  leaving  the 
penetrator  to  intercept  the  target. 
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Figure  1-3  PHALANX  MK  149  Round 
The  key  characteristics  of  the  gun  are  defined  in  Table  1. 


PARAMETER 

CHARACTERISTIC 

GUN 

Fire  Rate 

3,000  or  4,500  rounds/minute 

Ammunition  Magazine  Capacity 

1,556  rounds  (including  chuting) 

Reload  Time 

<  1 5  minutes 

PAC  Burst  Length 

Continuous,  1.2  seconds  or  2.0  seconds 

ROUND  (MK  149) 

Muzzle  Velocity 

>  3600  feet/second  (1027  m/s) 

Velocity  Decrease 

0.2  feet/second  per  foot  (0.06  m/s  per  m) 

Round  Weight 

0.6  pounds  (0.27  kg) 

Penetrator  Weight 

>  1000  grains 

Caliber 

20  mm 

Table  1-1,  Key  Characteristics  of  the  PHALANX  Gun  System 
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B.  MOTIVATION 

In  today's  littoral  operating  environment  the  threat  is  unclear  and  unpredictable, 
and  Navy  ships  are  in  constant  threat  from  anti-surface  missiles  (ASM's).  These  missiles 
are  getting  faster,  smaller,  smarter  and  are  available  to  almost  any  country  in  the  world. 
The  PHALANX  close-in  weapon  system  is  a  very  important  part  of  the  anti-ship  missile 
defense  structure  of  the  U.S.  Navy.  The  more  that  is  understood  about  parameters  that 
control  PHALANX  gun  dispersion  performance,  the  more  the  Navy  can  increase  its 
performance  and  therefore  ship  survivability. 

The  data  collected  on  a  firing  PHALANX  gun  at  the  Naval  Air  Weapon  Station, 
China  Lake  CA,  will  help  in  understanding  the  relationship  between  the  barrel  tip 
movement  and  dispersion.  The  objective  of  this  research  is  to  correlate  accelerometer 
data  and  dispersion  data.  The  ultimate  goal  is  to  apply  the  results  of  this  study  to 
PHALANX  systems  in  the  fleet  and  determine  their  dispersion  value. 

This  research  was  conducted  in  support  of  the  Naval  Surface  Warfare  Center 
(NSWC),  Port  Hueneme  Division,  testing  and  evaluation  of  the  PHALANX  close-in 
weapon  system. 
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II.  DESCRIPTION  OF  FIELD  EXPERIMENTS 


A.  INTRODUCTION 

The  Naval  Air  Weapon  Station  (NAWS)  at  China  Lake,  CA,  is  one  of  the  Navy's 
weapons  testing  facilities.  It  is  used  for  testing  and  evaluation  of  new  and  current  weapon 
systems.  It  is  at  NAWS  that  the  Naval  Surface  Weapon  Center  (NSWC)  Port  Hueneme 
Division  conducts  all  testing  of  the  PHALANX  close-in  weapon  system.  NSWC 
conducts  many  performance  evaluations  of  PHALANX.  They  include  performance 
versus  live  targets,  sonic  and  subsonic,  and  tracking  exercises  versus  numerous  targets. 
The  facility  is  also  equipped  with  round  dispersion  analysis  equipment. 

1.  Equipment  setup 

Data  for  dispersion  analysis  were  collected  on  the  PHALANX  CIWS  system  at 
the  NAWS.  (Figure  2-1)  The  CIWS  system  used  was  a  Block  1  Baseline  2  with  HOL 
(High  Order  Language).  It  was  part  of  a  CTE  (Contractor  Testing  &  Evaluation) 
evaluation  phase  of  hardware  and  software  testing.  During  a  CTE,  the  contractor, 

Hughes  Corporation,  is  allowed  to  use  government  facilities  for  research  and 
development. 

Two  trips  were  made  to  China  Lake,  the  first  one  on  9  May  1994  and  the  second 
one  on  31  May  1994.  Eleven  PAC  (Pre-Action  Calibration)  fires  were  conducted  on  each 
day  at  low  and  high  firing  rates,  firing  between  100  and  170  rounds,  respectively.  On 
each  firing,  accelerometer  data  were  collected  at  various  positions  on  the  PHALANX  gun 
system.  PACs  were  conducted  with  the  prototype  muzzle  restraint,  production  muzzle 
restraint,  and  with  no  muzzle  restraint.  The  accelerometer  data  were  used  to  determine 
the  gun  movement  during  a  firing.  Dispersion  data  were  also  collected  but  only  on  the 
second  trip  to  China  Lake.  The  dispersion  data  were  taken  to  correlate  against  the  gun 
movement,  specifically  lateral  barrel  tip  deflection  with  the  production  muzzle  restraint 
installed. 
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Figure  2-1  PHALANX  Gun  System  at  China  Lake,  CA 
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a.  Dynamic  Signal  Analyzer  Model  SA  390 

In  order  to  conduct  a  study  and  understand  the  movement  of  the 
PHALANX  gun  system,  accelerometer  data  were  collected  using  the  Scientific  Atlanta 
SA  390  Dynamic  Signal  Analyzer.  [Ref  3]  The  SA  390  was  donated  to  the  Naval 
Postgraduate  School  by  Scientific-Atlanta.  The  SA  390  is  a  full-featured  8  channel,  0- 
100  kHz  signal  analyzer  featuring  16-bit  AID  conversion  with  32-bit  floating  point 
processing  for  high  speed  digital  signal  processing.  The  SA  390  provides  real-time 
processing  of  data  for  which  acquisition,  processing  and  display  can  take  place  with  no 
loss  of  data.  It  has  selectable  resolution  of  200, 400, 800  or  1600  spectral  lines.  The  SA 
390  was  the  signal  analyzer  of  choice  because  of  its  eight  channels  and  extended  recorder 
capability. 

b.  PCB  Quartz  Shear  Mode  ICP  Accelerometers 

The  accelerometer  data  were  collected  using  PCB  Quartz  Shear  Mode  ICP 
Accelerometer  Series  353.  [Ref  4]  The  Series  353  accelerometers  employ  built-in,  solid 
state  microelectronics  for  conditioning  the  signal  generated  by  the  naturally  piezoelectric, 
quartz  crystal,  sensing  element.  This  type  of  sensor  is  designated  ICP  which  stands  for 
Integrated  Circuit  Piezoelectric.  The  built-in  ICP  microelectronics  incorporates  a 
MOSFET  transistor  that  requires  excitation  power.  This  is  supplied  by  a  regulated 
constant  current  source  of  typically  4  milliamps. 

In  order  to  conduct  a  complete  study  of  motion,  sets  of  accelerometers  were 
placed  on  three  key  locations  on  the  gun;  the  muzzle  restraint,  the  mid-house  bearing  and 
the  aft  bearing.  They  were  mounted  on  tri-axial  cubes.  The  muzzle  restraint  and  aft 
bearing  accelerometers  were  mounted  with  their  sensitive  axis  oriented  horizontally  and 
vertically.  The  mid-house  accelerometers  were  mounted  along  the  horizontal,  vertical 
and  axial  axes.  Figures  2-2, 2-3, 2-4  show  these  mounts. 
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Figure  2-4  Aft  Bearing  Accelerometers 
Table  2-1  summarizes  the  accelerometer  set  up  for  9  May  1994. 


Table  2-1,9  May  1994  Accelerometer  Set  Up 
c.  Oehler  System  82 

The  facility  at  NAWS  is  equipped  with  an  Oehler  System  82.  The  Oehler 
system  is  used  to  collect  dispersion  data  on  a  firing  gun.  It  consists  of  a  frame  with  a 
microphone  at  each  comer.  For  a  dispersion  firing,  the  PHALANX  gun  is  aimed  at  the 
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center  of  the  frame,  thus  forming  an  apparent  target  for  the  PHALANX  supersonic 
projectiles.  The  Oehler  system  detects  where  within  the  frame  the  bullet  passed,  and  with 
the  machine-gun  chronograph,  it  records  bullet  position  and  time  data.  It  calculates  the 
rate-of-fire  for  each  round.  The  Oehler  System  82  merges  measurement,  recording  and 
analytical  functions  required  in  ballistic  testing  into  a  single  integrated  system.  Oehler 
dispersion  data  will  be  used  to  correlate  the  dispersion  against  the  gun  movement  data 
obtained  from  the  accelerometers.  Figure  2-5  shows  the  Oehler  frame  and  battenboard 
system  at  NAWS  China  Lake.  The  frame  is  six  feet  in  diameter  and  the  system  accuracy 
is  typically  0.1  inches.  [Ref  5]  The  battenboard  provides  an  immediate  picture  of  the 
dispersion.  Appendix  A  contains  the  output  of  the  Oehler  System  for  the  firing  three. 
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Figure  2-5  Oehler  Frame  and  Battenboard  System  at  China  Lake,  CA 
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B.  EXPERIMENTAL  SETUP  FOR  9  MAY  1994 


On  the  first  trip  to  China  Lake,  accelerometer  data  were  collected  on  a  firing  gun 
with  the  prototype  muzzle  restraint  (Figure  2-6),  the  production  muzzle  restraint  (Figure 
2-7),  and  w'ith  no  muzzle  restraint.  The  muzzle  restraints  were  designed  by  Hughes  as  an 
attempt  to  control  dispersion.  The  production  muzzle  restraint  is  the  restraint  approved  as 
a  system  modification  for  the  PHALANX  system.  It  is  being  produced  and  provided  to 
the  U.S.  Navy  as  an  ordnance  modification  for  its  close-in-weapon  systems. 


Figure  2-6  Prototype  Muzzle  Restraint 


Figure  2-7  Production  Muzzle  Restraint 

The  9  May  trip  to  China  Lake  was  used  as  a  learning  experience  for  data 
collection.  Terrible  ground  loop  problems  were  experienced  while  collecting  data  on  the 
PHALANX  gun  system.  This  occurred  because  the  two  accelerometers  used  at  the 
muzzle  restraint  were  not  ground  isolated.  For  later  data  collection,  this  problem  was 
corrected  by  using  J-series  (ground  isolated)  accelerometers.  These  accelerometers 
prevent  ground  loops  by  isolating  the  accelerometer  case  from  the  base. 

For  data  collection,  the  intent  was  to  use  the  firing  signal  from  PHALANX  gun 
system  as  a  trigger  for  the  SA  390 ,  but  as  seen  in  Figure  2-8,  the  signal  was  very  erratic 
and  a  quick  solution  was  required.  (The  erratic  firing  signal  was  identified  by  the 
contractor  as  due  to  modifications  in  the  new  HOL  system  software.)  A  manual  jury 
rigged  trigger  was  used.  The  SA  390  triggers  on  a  TTL  high  to  low  transition.  A  one 
kilohertz  square  wave  was  generated  at  TTL  levels  using  a  Philips  Model  PM  3384  DSO. 
[Ref  6]  It  was  manually  applied  about  half  a  second  before  firing.  A  button  trigger  was 
built  for  the  second  trip.  This  trigger  was  a  voltage  divider  with  a  9  volt  battery.  The 
voltage  divider  applied  a  TTL  high  to  low  transition  w'hen  the  button  was  pushed.  The 
button  was  pushed  when  there  w’as  one  second  left  in  the  countdown.  The  signal  analyzer 
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commenced  collecting  data  immediately  prior  to  the  gun  firing.  This  was  done  in  order 
to  see  and  study  the  beginning  of  the  gun  firing. 


A-CompTime-L1 

Firing  Signal  Fire  3  (4500)  9  May  1994 
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Figure  2-8  Firing  Signal  for  the  Third  Firing  9  May  1994 


The  SA  390  sampling  rate  and  time  record  can  be  adjusted  for  any  need.  The 
sampling  rate  is  determined  by  the  selected  frequency  span  (4  kHz)  using  the  equation 

2.56  X  FrequencySpan  =  SamplingRate 
2.56  X  4kHz  =  10,240  samples / channel f sec 

For  the  first  trip  to  China  Lake,  it  was  predetermined  that  3.2  seconds  would  be 
sufficient  recording  time.  The  SA  390's  extended  recorder  time  is  determined  by  the 
sampling  rate  and  memory  allocated.  With  eight  channels,  the  number  of  samples  per 
channel  was  30,768  and  the  total  memory  allocated  per  firing  was  256k.  Figure  2-9 
shows  that  3.2  seconds  was  not  sufficient  to  collect  the  ramp-down  after  firing.  For  the 
second  trip  the  time  records  were  increased  to  four  seconds,  which  proved  to  be  sufficient 
to  record  a  complete  PAC  fire.  Increasing  the  extended  recorder  time  increased  the 
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number  of  samples  per  channel  to  40,960.  No  dispersion  data  were  collected  on  9  May, 
therefore  no  correlation  between  the  accelerometer  data  and  dispersion  data  could  be 
completed. 


D“CompTimc-L4 
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Figure  2-9  Muzzle  Restraint  Acceleration  (Horizontal)  for  9  May  1994 

On  9  May,  the  PHALANX  magazine  was  fully  loaded  with  1500  rounds.  This 
allowed  for  eleven  PAC-fires.  The  firing  rate  was  varied  between  3000  RPM  and  4500 
RPM.  The  amount  of  rounds  fired  varied  between  100  and  170.  Data  were  collected 
with  the  prototype  muzzle  restraint,  the  production  muzzle  restraint  and  without  muzzle 
restraint.  Tables  2-2, 2-3  and  2-4  summarizes  the  firing  runs  conducted  on  9  May  1994. 


Firing  Run  Number 

Firing  Rate  RPM 

Rounds  Fired 

1 

3000 

100 

2 

4500 

170 

3 

4500 

170 

Table  2-2,  9  May  Firing  runs  made  with  prototype  muzzle  restraint 
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Firing  Run  Number 

Firing  Rate  RPM 

Rounds  Fired 

4 

3000 

100 

5 

4500 

160 

6 

3000 

120 

7 

4500 

160 

8 

3000 

110 

9 

4500 

No  Data  Available 

Table  2-3,  9  May  Firing  rans  made  with  production  muzzle  restraint 


Firing  Run  Number 

Firing  Rate  RPM 

Rounds  Fired 

10 

3000 

No  Data  Available 

11 

4500 

No  Data  Available 

Table  2-4,  9  May  Firing  runs  made  with  no  muzzle  restraint 

From  the  experience  learned  on  9  May,  it  was  decided  that  a  second  trip  to  China 
Lake  would  be  extremely  beneficial.  This  trip  was  conducted  on  31  May  1994.  This  time 
allowed  for  the  ordering  of  two  additional  J-series  accelerometers  and  preparing  better 
mounts  for  the  tri-axial  cubes. 

C.  EXPERIMENTAL  SETUP  FOR  31  MAY  1994 

The  same  experimental  set  up  was  used  for  the  second  trip  to  China  Lake, 
although  accelerometer  data  were  only  obtained  from  the  production  muzzle  restraint. 

The  following  changes  were  implemented  from  the  first  trip:  bullet  dispersion  data  were 
collected  with  the  Oehler  frame,  a  push-button  trigger  was  used,  and  ground-isolated 
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accelerometers  were  used  for  all  seven  locations  on  the  gun.  Also,  temperature 
measurements  were  made  at  the  muzzle  restraint  before  and  after  each  firing  to  see  if 
there  was  any  variation  in  the  dispersion  with  temperature.  This  study  does  not  take  into 
account  the  temperature  effects. 

The  SA  390  was  configured  to  record  four  seconds  of  data.  Figure  2-10,  shows 
that  this  was  sufficient  time  to  record  the  ramp-up  and  ramp-down  of  the  gun  firing. 
Again,  eleven  PAC  fires  were  conducted.  Table  2-5  provides  a  summary.  The 
temperatures  listed  for  each  firing  were  the  measured  temperature  before  and  after  each 
firing.  The  data  collected  on  31  May  1994  were  analyzed  for  research. 
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Figure  2-10  Muzzle  Restraint  Acceleration  (Horizontal)  for  31  May  1994 
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Firing  Run 

Firing  Rate 

Rounds  Fired 

Dispersion 

Temperature  (C) 

Before  Firing  After  Firing 

1 

4500 

170 

0.93 

38.8 

42.8 

2 

3000 

110 

0.89 

44.3 

48.8 

3 

4500 

170 

0.93 

50.1 

55.0 

4 

3000 

110 

1.06 

56.1 

58.2 

5 

3000 

120 

1.06 

59.6 

61.9 

6 

4500 

170 

1.18 

61.3 

61.8 

7 

4500 

160 

0.91 

62.5 

64.3 

8 

4500 

170 

0.89 

63.9 

66.2 

9 

4500 

170 

1.14 

66.2 

68.7 

10 

4500 

160 

1.06 

68.0 

68.9 

11 

3000 

50 

N/A 

57.3 

58.0 

Table  2-5,  Firing  Runs  for  31  May  1994  (all  made  using  the  production  muzzle  restraint) 

Figures  2-11  through  2-18  show  the  output  for  all  eight  channels  of  the  SA  390  for  firing 
run  three.  These  data  were  used  in  all  subsequent  analyses. 
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B-CompTime-L2 

Muzzle  Restraint  Acceleration  (Horizontal)  Fire  3 


Figure  2-11  SA  390  Trace  for  the  Muzzle  Restraint  Acceleration  (Horizontal) 


C-CompTrme-L3 

Muzzle  Restraint  Acceleration  (Vertical) 


Figure  2-12  SA  390  Trace  for  the  Muzzle  Restraint  Acceleration  (Vertical) 
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Figure  2-13  SA  390  Trace  for  the  Mid-House  Bearing  Acceleration  (Horizontal) 


Figure  2-14  SA  390  Trace  for  the  Mid-House  Bearing  Acceleration  (Vertical) 
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F>CompTime*L6 

Mid-House  Bearing  Acceleration  (Axial) 


Figure  2-15  SA  390  Trace  for  the  Mid-House  Bearing  Acceleration  (Axial) 
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Figure  2-16  Aft  Bearing  Acceleration  (Horizontal) 
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H-CompTfme-L8 

Afl  Bearing  Acceleration  (Vertical) 


Figure  2-17  SA  390  Trace  for  the  Aft  Bearing  Acceleration  (Vertical) 


A-C.ompTIme-LI 

Firing  Signal  Fire  3  (4500)  31  May  1994 


Figure  2-18  Firing  Pin  Signal 
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III.  DISPERSION  DATA  ANALYSIS 


Dispersion  is  directly  related  to  the  hit-probability  of  any  gun/ammunition  system. 
Many  factors  affect  dispersion,  but  the  most  important  are  the  wear  on  the  system,  the 
type  of  round  fired,  and  the  firing  rate.  The  firing  rate  has  a  major  impact  on  dispersion. 
For  example,  when  firing  at  200  rounds/minute  or  higher  the  dispersion  increases  by  50 
per  cent  above  the  single  shot  figure.  [Ref  7]  Other  factors  that  affect  dispersion  are 
wind,  tracking  errors  between  the  sight  and  the  gun,  trunnion  tilt,  and  recoil  forces.  [Ref 
8]  When  a  gun  is  manufactured,  it  will  be  tested  for  accuracy  and  dispersion  against  a 
test  specification.  A  worst-case  weapon,  which  only  just  passes  its  test  specification  as 
new,  would  have  nearly  twice  the  dispersion  of  a  typical  production  gun.  Also,  since  all 
guns  wear  during  their  life,  the  dispersion  will  increase  by  up  to  20  per  cent  before 
refurbishment  is  considered  necessary.  [Ref  9] 

The  M61A1  gun  was  initially  designed  for  aircraft  use.  When  used  on  aircraft, 
high  dispersion  is  effective.  It  gives  the  aircraft  a  higher  probability  of  a  kill  against 
another  fast  moving  target.  When  the  gun  was  adopted  for  the  Navy's  close-in  weapons 
system,  it  performed  satisfactorily  against  slow  sea-skimming  anti-ship  missiles.  The 
radar  technology  was  still  under  development,  and  dispersion  had  not  been  identified  as  a 
problem.  As  the  radar  improved,  it  began  to  monitor  erratic  dispersion  patterns.  These 
patterns  cause  the  CIWS  closed-loop  spotting  to  pull  the  projectile  stream  off  the  target. 
Radar  technology  has  improved  dramatically  and  surpassed  the  performance  of  the  gun. 
Today  the  anti-ship  missiles  are  faster  and  smarter.  Therefore,  PHALANX,  the  last 
weapon  in  the  defense  in  depth  concept,  must  destroy  them  at  the  furthest  possible  range. 
The  high  dispersion  leaves  a  deficiency  in  a  ship's  defensive  capability. 
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A.  DISPERSION  DATA 

1.  Firing  Rate 

Dispersion  data  was  acquired  on  31  May  1994  for  ten  of  the  eleven  firings.  The 
set  up  was  controlled  by  Yuji  Wilson,  the  Project  Manager  and  our  sponsor  for  the 
PHALANX  gun  system.  The  dispersion  is  determined  when  the  firing  rate  is  steady.  It 
takes  an  average  of  ten  to  twelve  rounds  for  the  gun  to  ramp  up  to  the  desired  firing  rate. 
The  gun  sometimes  overshoots  the  desired  firing  rate,  but  it  settles  back  down.  This  was 
seen  when  the  firing  rate  was  decreased  from  4500  RPM  to  3000  RPM.  Figure  3-1 
shows  these  features  for  firing  run  number  two. 


Fire  2  3000  RPM  -  102  Rounds 


Figure  3-1  Firing  Rate  for  the  Second  Firing,  31  May  1994 


The  Oehler  82  System  records  the  time  of  each  round  arrival  at  the  frame  and  its 
location.  From  this,  it  determines  the  gun  firing  rate.  The  two  times,  0.45  s  and  1.79  s, 
correspond  to  when  the  gun  achieved  steady  state  and  when  it  commenced  to  ramp  down. 
Figure  3-2  shows  the  firing  rate  for  the  third  firing. 
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Fire  3  4500  RPM  - 161  Rounds 
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Figure  3-2  Firing  Rate  for  the  Third  Firing,  31  May  1994 

The  marks  at  around  4.5  seconds  are  flyers.  Flyers  can  be  sabot  petals  or 
projectile  debris  that  hits  one  of  the  microphones  directly  on,  or  a  bullet  that  hits  outside  a 
certain  radius.  Flyers  are  disregarded  when  determining  dispersion.  In  order  to  estimate 
a  correct  value  for  dispersion,  the  rounds  with  a  firing  rate  10  percent  above  and  below 
the  preset  firing  rate,  are  eliminated.  These  rounds  account  for  the  ramp-up  and  ramp- 
down  of  the  gun.  If  there  are  flyers,  the  Oehler  computer  calculates  a  radial  value  for  sub 
bursts.  A  sub  burst  includes  rounds  at  steady  firing  up  to  the  first  flyer.  A  radial  value  is 
calculated  for  each  sub  burst.  The  process  continues  until  all  the  rounds  within  the  steady 
rate  are  accounted  for.  The  radial  value  for  each  sub  burst  is  multiplied  by  a  factor 
determined  by  the  number  of  rounds  in  the  sub  burst.  The  radial  values  are  summed  and 
divided  by  the  sum  of  the  factors.  This  value  is  multiplied  by  0.707  to  determine  the 
estimated  dispersion  value.  Table  3-1  summarizes  the  firing  runs  for  31  May  1994. 
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Firing  Run 

Firing  Rate  (RPM) 

Rounds  Fired 

Dispersion  (mrads) 

1 

4500 

156 

0.93 

2 

3000 

102 

0.89 

3 

4500 

159 

0.93 

4 

3000 

101 

1.06 

5 

3000 

101 

1.06 

6 

4500 

160 

1.18 

7 

4500 

157 

0.91 

8 

4500 

156 

0.89 

9 

4500 

154 

1.14 

10 

4500 

154 

1.06 

11 

3000 

50 

N/A 

Table  3-1  Firing  Runs  for  31  May  1994 

There  is  a  disparity  between  Table  2-5  and  Table  3-1  in  the  number  of  rounds 
fired.  In  Table  2-5,  the  number  of  rounds  fired  was  determined  by  the  firing  pin  signal. 
As  stated  earlier,  this  signal  was  extremely  erratic.  The  number  of  rounds  fired  listed  in 
Table  3-1  was  determined  by  the  Oehler  Frame,  and  are  the  correct  values. 

2.  Dispersion 

Figure  3-3  is  a  plot  of  the  dispersion  pattern  for  firing  run  3.  The  fact  that  the 
pattern  is  not  centered  at  (0,0)  is  not  significant.  The  gun  barrels  are  aligned  to  the  center 
of  the  Oehler  Frame,  which  is  eighty  three  feet  away.  It  is  difficult  to  get  the  barrels 
aimed  at  the  exact  center  of  the  frame.  This  is  easily  corrected  by  subtracting  the  mean. 
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Dispersion  Pattern  for  the  Third  Firing 
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Figure  3-3  Dispersion  Pattern  for  the  Third  Firing 

In  order  to  cross  correlate  the  displacement  of  the  muzzle  restraint  against  the 
dispersion,  the  dispersion  data  was  transferred  into  MATLAB.  The  array  of  dispersion 
data  had  to  be  expanded  from  157  elements  long  to  40960  in  order  to  match  the  size  of 
the  accelerometer  data.  The  first  column  of  the  dispersion  matrix  contains  the  best 
estimate  of  the  time  of  the  first  round  fired.  This  value  was  determined  using 
CROSSCORR.m  (Appendix  B).  The  procedure  is  detailed  in  the  next  chapter.  The 
second  and  third  column  contain  the  horizontal  and  vertical  dispersion  data. 
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rV.  ACCELEROMETER  DATA  ANALYSIS 


Due  to  time  constraints,  it  was  decided  to  concentrate  the  efforts  of  signal  analysis 
on  one  firing  run.  Firing  run  number  three  of  3 1  May  1994  was  chosen  for  no  particular 
reason. 

The  SA  390  Dynamic  Signal  Analyzer,  although  a  powerful  tool,  is  limited  in  its 
data  analysis  capabilities.  Therefore,  in  order  to  conduct  an  analysis,  it  was  necessary  to 
transfer  the  data  to  a  powerful  technical  computing  environment  for  high-performance 
numeric  computation  and  visualization  such  as  MATLAB.  [Ref  10] 

A.  DATA  CONVERSION 

The  data  stored  in  the  SA  390  is  in  a  raw  ATT  32C  DSP  binary  format.  [Ref  1 1] 

It  was  required  to  convert  the  data  to  ASCII  format  for  use  by  PC's.  The  S  A  390  does  not 
have  the  capability  of  exporting  extended  records.  Initially  there  was  no  choice  in  the 
method  of  converting  and  exporting  the  data.  The  only  approach  was  to  load  the 
extended  record,  and  save  each  successive  four  hundred  millisecond  time  trace  as  a  text 
file.  The  ten  text  files  were  then  reconstructed  in  a  text  editor  such  as  Brief.  Since  data 
were  collected  on  twenty  firings,  this  process  was  tedious  and  long.  With  the  help  of 
Robert  Mihata  of  Scientific  Atlanta,  an  executable  code,  XRC2MAT,  was  eventually 
obtained.  XRC2MAT  extracts  the  data  from  an  extended  record  and  saves  it  to  another 
file  in  ASCII  format.  When  this  file  is  loaded  in  MATLAB,  each  channel  of  data 
collected  represents  a  column  in  a  matrix.  This  resulted  in  a  40960  by  8  matrix. 

With  the  data  converted  and  imported  into  MATLAB,  the  analysis  of  the  data 
began.  The  goal  of  the  data  analysis  was  to  correlate  the  barrel  tip  lateral  displacement 
with  the  dispersion.  The  first  objective  was  to  properly  display  the  data  in  MATLAB  and 
convert  it  from  volts  to  meters  per  second^2.  With  the  help  of  Mike  Hatch,  the  code 
GTESTFFT.m,  (Appendix  C),  was  written.  GTESTFFT.m  originally  converted  the 
acceleration  from  volts  to  inches  per  second'^2  and  plots  the  acceleration  versus  time. 
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This  procedure  was  modified  to  convert  the  acceleration  from  volts  to  meters  per 
second''2.  It  creates  the  time  scale  by  dividing  the  number  of  samples  by  the  sampling 
rate  (10,240  samples/second).  Figure  4-1  shows  the  reconstructed  record  of  the  aft 
bearing  accelerometers.  Figure  4-2  shows  the  acceleration  for  the  muzzle  restraint 
accelerometers.  Figure  4-3  shows  the  acceleration  for  the  mid  housing  bearing 
accelerometers.  GTESTFFT.m  also  computes  and  plots  the  normalized  discrete-time 
Fourier  transform  (Figure  4-4)  by  using  the  following  equation  from  the  Signal 
Processing  Toolbox  in  MATLAB  [Ref  12,  p.  2-90  2-92] 

Pn^abs(ffiix))*lllength{x)  (1) 


Aft  Bearing  Acceleration  (Vertical) 


1 0000 1 - 1 - r 


Aft  Bearing  Acceleration  (Horizontal) 


15000 


Time,  sec  -  Fire  3 


Figure  4-1  Aft  Bearing  Acceleration 
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Muzzle  Restraint  Acceleration  (Horizontal) 


Figure  4-3  Mid  Housing  Bearing  Acceleration 
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Normalized  FFT  of  Muzzle  Restraint  Acceleration  (Horizontal) 
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Figure  4-4  Normalized  FFT  of  Muzzle  Restraint  Acceleration  (Horizontal) 


B.  DATA  ANALYSIS 

1.  Determining  the  Time  of  the  First  Round 

In  order  to  commence  the  correlation  between  measured  acceleration  and 
dispersion,  it  was  required  to  time-align  the  dispersion  data  and  the  accelerometer  data. 

If  the  firing  pin  signal  had  not  been  so  erratic,  then  it  would  have  been  used  for  this 
purpose,  and  time  aUgnment  would  have  been  a  trivial  problem.  Determining  the  time 
correspondence  between  the  acceleration  and  dispersion  data  records  was  completed  in 
two  steps.  First,  from  the  acceleration  plots  it  was  determined  that  the  first  round  fired 
occurred  within  the  first  12228  samples.  Using  THRESHOLD.m  (Appendix  D),  a  limit 
value  was  used  to  determined  the  first  occurrence  of  an  acceleration  magnitude  exceeding 
a  threshold  value.  THRESHOLD.m  returned  the  bin  number  when  the  limit  was 
exceeded.  From  acceleration  values  around  the  exceeded  acceleration  magnitude,  the 
threshold  value  was  adjusted  until  the  bin  number  corresponded  to  the  maximum 
acceleration  magnitude.  For  data  of  firing  run  three,  this  gave  a  time  after 
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commencement  of  accelerometers  data  acquisition  of  0.9879  seconds.  Due  to  the  process 
of  resetting  the  threshold  value,  this  time  was  considered  a  rough  estimate  of  the  first 
round  firing  time. 

In  order  to  obtain  the  best  estimate,  it  was  decided  to  cross  correlate  the  dispersion 
and  accelerometer  data.  For  this  step  the  XCORR.m  function  in  the  Signal  Processing 
Toolbox  [Ref  13,  p.  2-226  2-228]  was  used.  To  cross  correlate  using  XCORR.m  the  data 
files  had  to  be  the  same  size.  As  stated  earher,  the  dispersion  array  was  increased  to  a 
40960  element  array.  XCORR.m  uses  FFT's  to  cross  correlate,  and,  due  to  the  size  of  the 
files,  there  were  memory  problems.  A  new  cross  correlating  function  had  to  be  created. 
CROSSCORR.m  correlates  for  300  values  around  the  rough  estimated  time,  50  prior  and 
250  after.  It  uses  the  following  equation  to  cross  correlate  between  the  dispersion  data 
matrix  x(n,l)  and  the  sum  of  the  squares  of  the  horizontal  and  vertical  muzzle  restraint 
accelerations,  batacc. 


..  4096a-(m| 

xcor= -  Y)  xcor(7n,l)+x(M,l)*6atacc(n+m,l)  (2) 

40,960- 1 OT I  tr 


Figure  4-5  shows  the  cross  correlation  between  the  accelerations  and  dispersion. 
The  time  closest  to  the  rough  estimate  with  the  largest  correlation  was  used  as  the  best 
estimate  of  when  the  PHALANX  gun  began  to  fire.  This  value  was  0.9836  seconds.  It  is 
presumed  that  this  time  is  more  accurate. 
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Cross  Correlation  Between  RredataS  &  batacc 


Figure  4-5  Results  of  Cross  Correlation  between  acceleration  and  dispersion 

The  results  from  THRESHOLD.m  provided  some  insight  into  the  PHALANX  gun 
movement  during  a  firing.  It  was  determined  that  the  gun  begins  to  move  first  at  the  aft 
bearing.  This  is  logical,  since  at  the  start  of  the  firing  sequence,  a  round  is  loaded  into  the 
barrel.  As  the  loaded  barrel  rotates  to  the  firing  position,  twelve  to  fourteen  degrees  from 
top  dead  center,  it  induces  movement  at  the  barrel  tip.  The  mid  housing  remains 
relatively  motionless  until  the  firing  signal  is  applied  to  the  round  and  the  round 
commences  to  travel  down  the  barrel. 

Figure  4-6  shows  the  muzzle  restraint  acceleration  with  the  best  estimate  of  the 
first  firing  time. 
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Figure  4-6  Muzzle  Restraint  Acceleration  (Horizontal)  versus  the  Best  Estimated  First 
Round  Time 

2.  Barrel  Tip  Lateral  Displacement  from  Acceleration 

Once  the  time  values  for  the  accelerometer  data  were  established,  we  commenced 
to  find  the  muzzle  restraint  movement.  Two  methods  were  used  to  find  and  compare  the 
lateral  displacement  at  the  muzzle  restraint:  direct  time  integration,  and  integration  using 
the  inverse  FFT. 

a.  Displacement  from  Direct  Time  Integration 

The  trapezoidal  numerical  integration  method  was  used.  The  trapezoidal 
rule  for  estimating  the  value  of 

b 

f  fix)  dx  (3) 

a 

is  based  on  approximating  the  region  between  the  graph  of  the^x)  and  the  x-axis  with 
trapezoids  of  equal  width.  The  trapezoids  have  the  common  base  h=(b-a)/n  and  the  side 
of  the  trapezoid  runs  from  the  x-axis  up  (or  down)  to  the  curve.  The  area  of  the 
trapezoids  add  up  to 
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T  -  ^  ^  *  Ax,))  h  .  ... 

*fK))  h  *  \(AxJ  ^Axji)  h 

(4) 

"  /(*I>  /(^a)  "  -  "  /(*„-!> 

»|  (/■(V  *  2  /(x,>  2  Ax,)  *  ...  *  2  /(*  )) 

where yfX(,j=^aj  and  Ax„)=f(b).  The  trapezoid  rule  estimates  the  integral  from  a  to  b. 

This  approach  was  used  to  obtain  the  displacement  from  the  acceleration  in  the  time 
domain.  The  following  equation  describes  the  integration  procedure. 

a.  for  i  =  1:40960 

'’i.i  =  v.  +  ^  (a.  *  a.  ^^)  ^  t  for  i  =  1:(40960-1) 

^  w) 

.  1  =  ^  (''i  ^  v.  ,)  A  t  for  i  =  l:(40960-2) 

where  a,  are  the  accelerations  values,  v,  are  the  velocity  values  and  X;  are  the  displacement 
values.  is  the  sampling  rate,  1/10240  sec. 

Figure  4-7  shows  the  results  of  double  integrating  the  muzzle  restraint 
acceleration  data.  From  the  figure,  the  muzzle  restraint  appears  to  have  about  a  four 
meter  displacement  range.  Also  note  that  at  time  zero  and  before  the  muzzle  restraint  is 
moving.  This  is  nonsense.  The  PHALANX  gun  is  only  displaced  by  millimeters. 
Therefore  the  displacement  of  interest  is  too  small  to  even  be  seen  in  this  figure.  The 
displacement  shown  is  dominated  by  a  large,  low  frequency  drift.  Figure  4-8  shows  the 
velocity  obtained  by  singly  integrating  the  accelerometer  data. .  The  actual  velocity  is 
seen  in  this  figure  as  the  "noise"  added  to  the  large,  low  frequency  trend. 
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Muzzle  Restraint  Displacement  (Horizontal) 


Figure  4-7  Muzzle  Restraint  Displacement  (Horizontal)  from  Direct  Time  Integration 


Muzzle  Restraint  Velocity  (Horizontal)  versus  Time 


Figure  4-8  Muzzle  Restraint  Velocity  (Horizontal)  from  Direct  Time  Integration 

It  is  speculated  that  the  large,  low  frequency  drift  is  due  to  temperatures  effects. 
Under  its  high  firing  rate,  the  PHALANX  barrels  tend  to  increase  rapidly  in  temperature. 
When  using  an  air  gun  to  heat  and  cool  an  accelerometer  to  recreate  this  effect,  there  was 
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a  drift  of  approximately  2  millivolts  per  degree  per  second.  As  seen  in  the  calibration 
charts  in  Appendix  E,  the  accelerometers  have  a  very  low  sensitivity  to  temperature,  but 
are  sensitive  to  high  gradient  temperature  changes.  Other  sources  of  the  low  frequency 
drift  could  be  the  high  impact  environment  of  the  firing  gun.  or  possibly  wind  effects  or 
the  electronics  (accelerometers  and  signal  analyzer). 

It  was  necessary  to  remove  the  large,  low  frequency  drift.  The  first  attempt  was 
to  use  the  DETREND .m  function  in  the  Signal  Processing  Toolbox.  [Ref  14,  p.  2-77] 
Using  DETREND.m,  the  linear  trend  and  mean  were  removed  from  the  acceleration. 
The  linear  trend  was  also  removed  from  the  velocity  and  from  the  displacement.  Note 
that  the  displacement  no  longer  begins  at  zero,  and  that  now  the  muzzle  restraint  appears 
to  have  a  displacement  range  of  two  and  one-half  meters.  (Figure  4-9)  This  is  a  step  in 
the  right  direction,  but  is  not  sufficient. 


Muzzle  Restraint  Displacement  (Horizontal) 


Figure  4-9  Detrended  Muzzle  Restraint  Displacement  (Horizontal)  from  Direct  Time 
Integration 


43 


b.  Displacement  from  the  FFT  of  the  Acceleration 
The  displacement  can  be  obtained  from  the  Fast  Fourier  Transform  of  the 
acceleration.  Dividing  the  FFT  of  the  acceleration  by  the  minus  frequency  squared  gives 
the  FFT  of  the  displacement.  Applying  the  Inverse  FFT  function  to  the  FFT  of  the 
displacement  gives  the  displacement  in  the  time  domain.  Equation  (6)  gives  the 
expressions  for  acceleration,  velocity  and  displacement  in  terms  of  the  coefficients  of 
their  fourier  transforms. 

a(t)  -  E  a, 

0 

v(0  =  E  V,  (6) 

0 

m  -  E  li 

0 

Equation  (7)  shows  the  relationship  between  the  coefficients  of  the  displacement 
transform  to  the  coefficients  of  the  acceleration  transform. 

HO  ■  ^  ■  i: 

at  0 

ait)  I,  (0^  (7) 

at  0 

1 

ur 

In  applying  this  relation  to  obtain  the  fourier  coefficients  of  the  displacement,  the 
fourier  coefficient  corresponding  to  zero  frequency  was  set  to  zero.  Figure  4-10  shows 
the  FFT  of  the  displacement  for  the  muzzle  restraint  in  the  vertical  axis.  The  peak  of 
77.29  Hertz  in  the  FFT  of  the  displacement  is  the  average  firing  rate. 
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Figure  4-10  FFT  of  Muzzle  Restraint  Displacement  (Vertical) 

Figure  4-11  shows  the  muzzle  restraint  displacement  obtained  using  this 

method. 


Muzzle  Restraint  Displacement  (Vertical) 


Time,  sec  -  Fire  3 

Figure  4-11  Muzzle  Restraint  Displacement  (Vertical)  from  the  Inverse  FFT 


c.  Comparison  of  Two  Integration  Methods 
Figure  4-12  shows  the  displacement  curves  obtained  from  the 
accelerometer  data  by  the  methods  described  in  the  preceding  sections.  Clearly  there  is  a 
disparity  between  the  two  results.  Several  procedures  to  improve  the  agreement  between 
the  two  were  investigated.  First,  the  acceleration  was  detrended  prior  to  time  integration. 
This  was  not  sufficient,  however.  Further  the  mean  was  removed  from  the  velocity  and 
displacement  curves  for  the  direct  time  integration.  This  procedure  resulted  in  excellent 
agreement  between  the  two  integration  methods,  and  therefore  either  one  could  be  used 
for  further  evaluation.  Figure  4-13  shows  these  results. 

Muzzle  Restraint  Displacement  (Vertical) 
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Muzzle  Restraint  Displacement  (Vertical) 


Figure  4-13  Comparison  of  the  Two  Integration  Methods  with  Detrending. 
d.  Removing  the  Low  Frequency  Trend 

In  order  to  derive  the  '"true"  barrel  tip  displacement  from  the  acceleration 
data,  a  means  of  removing  the  large,  low  frequency  in  the  measurements  was  required. 
Initially  it  was  thought  that  removing  the  trend  from  a  small  section  of  data  would  prove 
most  effective. 

Figure  4-14  shows  the  detrended  displacement  from  1.216  seconds  to  1.413 
seconds.  The  displacement  has  a  range  of  0.03  inches  or  0.762  millimeters.  These 
values  are  much  closer  to  the  actual  movement  of  the  gun.  The  integration  of  the  entire 
record  does  not  show  the  true  displacement.  By  forcing  the  integration  on  a  small  section 
when  the  gun  was  at  steady  firing  rate,  the  displacement  can  be  obtained.  Three  methods 
were  investigated  in  order  to  obtain  the  "true"  displacement  for  the  entire  record:  time  and 
frequency  windowing,  global  polynomial  approximation  and  cubic  spline  smoothing. 
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Figure  4-14  Detrended  Displacement  from  1.216  to  1.413  Seconds 

(1)  Time  and  Frequency  Windowing.  The  possibility  of  removing 
the  trend  from  a  time  windowed  record  by  a  high  pass  filter  in  the  frequency  domain  was 
investigated.  Two  window  functions  were  created,  one  in  the  time  domain,  TIMWIN.m 
(Appendix  F)  and  one  in  the  frequency  domain,  FREQWIN.m  (Appendix  G).  They 
allowed  for  a  smooth  transition  from  zero  to  a  desired  value.  These  window  functions  are 
composed  of  a  half  cosine  and  half  sine.  Figure  4-15  is  an  example  of  a  normalized  half 
cosine  and  half  sine  window  functions.  This  function  was  applied  to  the  acceleration  data 
in  the  time  domain  to  effect  a  time  window  or  to  FFT  data  in  the  frequency  domain  to 
effect  a  high-pass  filter.  (In  MATLAB,  the  FFT  includes  negative  frequencies;  small 
negative  frequencies  occur  at  the  end  of  the  FFT  record).  [Ref  15] 
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Figure  4-15  Window  Function 
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Figure  4-16  shows  the  time  domain  window  function  applied  to  the  muzzle 
restraint  acceleration  for  1.5  to  2  seconds. 


Muzzle  Restraint  Acceleration  (Horizontal) 


Time,  sec 


Figure  4-16  Time- Windowed  Muzzle  Restraint  Acceleration  (Horizontal)  From  1.5  to  2 
Seconds 

In  Figure  4-17  the  frequency  domain  window  was  applied  to  the  FFT  of  the 
acceleration,  then  it  was  divided  by  omega  squared,  the  result  is  the  windowed  FFT  of  the 
displacement.  The  FREQWIN.m  function  could  not  be  used  above  three  hertz  since  the 
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first  mode  of  vibration  is  3.4  Hz  as  identified  in  the  Finite  Element  Model.  [Ref  16]  The 
use  of  these  functions,  although  helpful  in  concentrating  on  one  region,  could  not  be 
applied  to  the  entire  record.  Windowing  had  to  be  abandoned  because  in  order  to  remove 
the  low  frequency  drift,  a  lot  of  the  data  had  to  be  excluded.  It  was  not  the  right 
approach. 


Muzzle  Restraint  Displacement  (Horizontal) 


Figure  4-17  Frequency-Windowed  FFT  of  the  Muzzle  Restraint  Displacement 
(Horizontal) 


(2)  Global  Polynomial  Approximation.  The  second  method  used 
to  remove  the  low  frequency  drift  was  polynomial  approximation.  Polynomials  are  the 
approximating  functions  of  choice  when  a  smooth  function  is  to  be  approximated  locally. 
[Ref  17,  p.  1-3]  A  relatively  small  data  set  can  be  approximated  smoothly  but  as  the 
interval  gets  larger,  the  degree  of  the  approximating  polynomial  may  be  unacceptably 
large.  For  the  muzzle  restraint  acceleration ,  a  twelfth  order  polynomial  did  not  provide 
a  good  enough  approximation.  It  was  not  possible  to  fit  polynomials  of  higher  order  to 
the  displacement  curve  due  to  memory  problems.  Figure  4-18  shows  a  twelfth  order 
polynomial  fit  of  the  displacement  curve  and  the  difference  between  the  two  on  the  same 
scale. 
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Mu27:te  Restraint  Displacement  (Horizontal) 


Figure  4-18  Twelfth  Order  Polynomial  Fit  to  the  Displacement  Data 


Figure  4-19  is  the  difference  from  the  fitted  polynomial  and  the  displacement 
curve.  A  twelfth  order  polynomial  is  not  a  good  approximation,  since  the  barrel  tip  does 
not  move  over  a  one  centimeter  range. 


Muzzle  Restraint  Displacement  (Horizontal) 


Figure  4-19  Barrel  Tip  Displacement  Extracted  from  Twelfth  Order  Polynomial  Fit  to  the 
Displacement  Data 


51 


(3)  Cubic  Spline  Approximation.  An  alternative  to  a  polynomial 
approximation  is  to  subdivide  the  interval  of  approximation  into  sufficiently  small 
intervals.  On  each  such  interval,  a  polynomial  of  relatively  low  degree  can  provide  a 
good  approximation.  A  cubic  spline  is  a  smooth  piecewise  polynomial,  where  the 
derivatives  of  each  polynomial  pieces  blend  smoothly.  I.J.  Schoenberg  coined  this  term 
since  a  twice  continuously  differentiable  cubic  spline  (with  a  sufficiently  small  first 
derivative)  approximates  the  shape  of  a  draftsman's  spline.  [Ref  18,  p.  1-3] 

The  CSAPS.m  command  in  the  Spline  Toolbox  [Ref  19,  p.  2-13], 
constructs  a  cubic  smoothing  spline  s  to  a  given  data  x,  y.  It  is  constructed  for  the 
specified  smoothing  parameter  p.  For  p=0,  this  is  a  least-squares  straight  line  fit  to  the 
data,  while  on  the  other  extreme,  i.e.,  for  p=l,  this  the  "natural"  cubic  spline  interpolant. 
[Ref  20,  p.2-13]  As  p  moves  from  0  to  1,  the  smoothing  spline  changes  from  one 
extreme  to  the  other.  The  interesting  range  of  p  is  often  near  l/(l+hVl6),  with  h  the 
average  spacing  of  the  data  abscissae.  For  the  accelerometer  data  h=l/10240,  or  the 
sampling  rate.  Since  the  data  are  uniformly  spaced,  the  smoothing  spline  is  expected  to 
follow  the  data  closely  for  p=l/(l+hVl6*10)  and  to  provide  satisfactory  smoothing  when 

p=l/(l-hhVl6/10). 

CSAPSXY.m,  Appendix  H,  is  code  that  applies  the  CSAPS.m 
command  to  the  muzzle  restraint  displacement  curves  for  the  horizontal  and  vertical 
directions.  It  plots  the  displacement  curve  obtained  from  direct  integration  along  with  the 
cubic  spline  approximation.  Figure  4-20  shows  the  results  with  a  smoothing  parameter 
p=l/(l+hVl6/10)=0.999....  The  two  curves  are  on  top  of  each  other. 
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Muzzle  Restraint  Displacement  (Horizontal) 


Figure  4-20  Cubic  Spline  Smoothing  of  the  Displacement  Data 

Figure  4-21  shows  a  small  section  of  the  curve  with  an  expanded  vertical  scale, 
and  it  shows  how  this  cubic  approximation  is  not  smooth  enough.  Figure  4-22  shows  the 
difference  between  the  displacement  curve  and  the  cubic  smoothing.  It  shows  a 
displacement  range  of  four  microns.  Since  the  barrel  tip  is  displaced  about  a  millimeter, 
the  cubic  smoothing  parameter  was  too  large,  and  the  data  followed  to  closely.  The 
smoothing  parameter  needed  to  be  adjusted. 
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Muzzle  Restraint  Displacement  (Horizontal) 


Figure  4-21  Zoom  of  the  Cubic  Smoothing  Approximation.  Obtained  using  CSAPS.m 
The  Smoothing  Parameter  p=0.999999999999418 


Time,  sec  -  Fire  3 


Figure  4-22  Muzzle  Restraint  Displacement  (Horizontal)  Extracted  Using  Cubic 
Smoothing  of  the  Displacement  Data 
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Through  trial  and  error  the  smoothing  parameter  was  decreased  in  order  to  achieve 
a  "more"  smooth  approximation.  The  smoothing  parameter  needed  to  be  decreased  by  a 
factor  of  one  billion,  for  an  appropriate  smoothing.  The  final  smoothing  parameter 
chosen  was  equal  to  0.999941795727021 1 .  Figure  4-23  shows  a  section  of  the  cubic 
smoothing  spline  approximation  on  the  displacement  curve. 


Muzzle  Restraint  Displacement  (Horizontal) 


Figure  4-23  Zoom  of  the  Cubic  Smoothing  Approximation.  Obtained  using  CSAPS.m. 
The  smoothing  parameter  p=0.9999417957207211. 

Figures  4-24  and  4-25  show  the  difference  between  the  cubic  smoothing 
approximation  and  the  displacement  curves  for  the  horizontal  and  vertical  direction  at  the 
muzzle  restraint.  The  range  of  displacement  obtained  is  very  reasonable.  It  is  expected 
for  the  vertical  direction  to  have  a  larger  range  in  the  displacement  than  in  the  horizontal 
direction 
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Time,  sec  -  Fire  3 


Figure  4-24  Muzzle  Restraint  Displacement  (Horizontal)  Extracted  Using  Cubic 
Smoothing  on  the  Displacement  Data 


Time,  sec  -  Fire  3 


Figure  4-25  Muzzle  Restraint  Displacement  (Vertical)  Extracted  Using  Cubic  Smoothing 
on  the  Displacement  Data 
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Examination  of  a  small  section  of  the  time  line  reveals  some  details  about  the 
PHALANX  while  firing.  Figure  4-26  is  a  plot  of  the  horizontal  displacement  for  one 
tenth  of  a  second.  It  shows  each  round  firing  and  the  subsequent  ringing.  Figure  4-27 
shows  the  vertical  displacement  for  the  same  period  of  time.  The  shots  are  not  as  clear  as 
in  the  horizontal  direction,  but  the  curve  gives  a  rough  indication  of  the  number  of  rounds 
fired  in  one  tenth  of  second  or  an  average  firing  rate  of  77  rounds  per  minute. 


Muzzle  Restraint  Displacement  (Horizontal) 


Figure  4-26  Muzzle  Restraint  Displacement  (Horizontal)  From  1.5  to  2  Seconds 


Figure  4-27  Muzzle  Restraint  Displacement  (Vertical)  From  1.5  to  2  Seconds 
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Figure  4-28  gives  an  insight  on  the  movement  of  the  gun  while  its  firing.  Besides 
the  displacement  being  larger  in  the  vertical  direction,  it  leads  the  horizontal  movement. 
This  can  be  explained  by  the  firing  action  of  PHALANX.  The  firing  sequence  begins 
about  twelve  to  fourteen  degrees  before  top  dead  center  creating  an  imbalance,  thus  the 
smaller  horizontal  displacement. 


Figure  4-28  Difference  Between  the  Horizontal  and  Vertical  Displacement 

C.  CORRELATION  OF  DISPERSION  DATA  AND  BARREL  TIP  LATERAL 
DISPLACEMENT 

Figure  4-29  shows  the  time  of  each  round  fired  obtained  from  the  dispersion  data 
shown  as  vertical  lines.  Also  shown  is  the  derived  vertical  barrel  tip  displacement.  It  is 
evident  from  the  figure  that  there  is  good  correspondence  between  firing  time  and 
maximum  vertical  displacement.  Figure  4-30  shows  the  time  of  each  round 
superimposed  over  the  horizontal  barrel  tip  displacement.  Note  the  lag  in  the  maximum 
displacement  after  each  firing.  From  these  figures  it  is  evident  that  when  a  round  is  fired 
most  of  the  force  is  in  the  vertical  direction.  But  due  to  the  construction  of  the  gun 
system,  horizontal  vibration  is  later  excited.  Also  note  the  best  estimated  time  of  bullet 
firing  obtained  from  CROSSCORR.m  does  not  quite  coincide  with  the  time  of  the 
maximum  vertical  displacement.  For  future  work  it  would  be  necessary  to  recalculate  the 
time  of  the  first  round  using  the  vertical  displacement  instead  of  the  acceleration. 
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Figure  4-29  Correlation  of  Barrel  Tip  Displacement  and  Dispersion  Data  (Vertical) 
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^q-4  Correlation  of  Displacement  and  Dispersion  Data  (Horizontal) 


Figure  4-30  Correlation  of  Barrel  Tip  Displacement  and  Dispersion  (Horizontal) 

Figures  4-31  and  4-32  were  made  to  investigate  further  the  correlation  of  barrel  tip 
motion  with  bullet  dispersion.  In  figure  4-31  is  plotted  the  vertical  dispersion,  in  radians, 
against  the  vertical  barrel  tip  displacement  occurring  at  the  estimated  time  of  buUet  firing, 
also  in  radians.  Figure  4-32  shows  the  same  for  the  horizontal  direction.  If  there  were 
good  correlation  between  the  two  displacements  plotted,  one  would  expect  the  data  points 
to  lie  close  to  a  line  of  slope  one.  There  is  no  such  correlation  agreement  in  Figures  4-31 
and  4-32.  This  suggests  that  further  refinement  of  the  procedure  to  extract  the  barrel  tip 
displacement  from  the  accelerometer  data  is  needed.  This  will  be  the  subject  to  follow-on 
investigation. 
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Figure  4-31  Dispersion  Correlation  (Vertical) 


Displacement,  Radians 


Figure  4-32  Dispersion  Correlation  (Horizontal) 
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V.  CONCLUSIONS 


The  objective  of  this  body  of  work  was  to  correlate  the  barrel  tip  transverse 
displacement  with  bullet  dispersion.  The  displacement  was  derived  from  accelerometer 
data  obtained  from  a  firing  PHALANX  gun  system  with  the  production  muzzle  restraint 
installed.  From  the  data  collection  and  analysis,  many  difficulties  were  uncovered  in 
measuring  the  motion  of  the  gun  and  its  influence  on  dispersion.  In  order  to  correlate 
dispersion  with  displacement,  it  was  necessary  to  time-align  the  measured  dispersion  with 
the  measured  acceleration.  This  study  also  gave  an  insight  into  the  PHALANX  gun 
motion  while  firing. 

Due  to  a  large,  low  frequency  drift  in  displacement  caused  by  a  combination  of 
electronics,  environment  and  temperature,  the  "true"  displacement  of  the  barrel  tip  was 
best  estimated  by  subtracting  out  a  smooth  trend  using  cubic  splines.  Although  there  was 
a  good  correlation  between  the  time  of  firing  obtained  from  the  dispersion  and  the  time  of 
maximum  vertical  displacement,  correlation  between  dispersion  data  and  extracted  lateral 
barrel  tip  displacement  was  disappointing. 

While  the  results  of  the  present  investigation  did  not  show  the  correlation  that  was 
hoped  for,  it  is  still  expected  that  good  correlation  can  be  found  between  bullet  dispersion 
and  properly  extracted  lateral  barrel  tip  displacement.  To  that  end,  the  following 
suggestions  are  made  for  further  investigation.  One,  a  potentially  significant  source  of 
error  is  the  correspondence  between  the  times  of  bullet  firings  and  accelerometer  time 
data.  What  is  needed  most  is  a  reliable  firing  pin  signal.  It  is  suggested  that  repairing  this 
signal  be  given  high  priority.  Two,  further  investigation  of  the  best  way  to  subtract  the 
large,  low  frequency  trend  in  the  displacement  obtained  from  the  accelerometer  data  is 
needed.  Cubic  spline  smoothing  appears  to  be  the  best  method,  as  was  demonstrated  in 
the  present  investigation.  However,  more  investigation  is  needed  of  how  to  optimally 
choose  the  smoothing  parameter.  Third,  if  some  means  could  be  found  to  measure  lateral 
barrel  tip  motion  using  non-contacting  transducers,  such  as  a  proximity  probe  or  a  laser 
vibrometer,  then  the  problem  of  the  large  low  frequency  drift,  if  it  is  related  to 
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temperature,  might  be  lessened. 


APPENDIX  A.  DISPERSION  DATA 


This  appendix  contains  the  dispersion  data  obtained  at  Naval  Air  Weapons 
Station,  China  Lake,  CA.  The  Oehler  System  82,  was  used  to  collect  the  dispersion  data 
from  a  firing  PHALANX  gun  system. 
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"OEHLER  SYSTEM  82,  VERSION  4.3,  PHONE  512-327-6900" 
DATE  =  94-05-31  ll:45HO.K.  ?_ 

ENTER  ID  (LESS  THAN  80  CHARACTERS) 

"053194d3, 4500spm,  lOOrds,  prod  muz  rest,  h7-hl2_" 

O.K.?_ 

SLOW  FIRE  ?N 
MACHINE  GUN  ?_ 

DISTANCE  =  10  ?  _ 

MINIMUM  VELOaTY  =  2000  ?  3300_ 

MINIMUM  VELOOTY  =  3300  ?  _ 

MAX  RATE  OF  FIRE  =  6000  ?  5000_ 

MAX  RATE  OF  FIRE  =  5000  ?  _ 

TARGET  SIZE  =  72.6  IN  ?  _ 

TEMPERATURE  =  70  DEG  F  ?  86_ 

TEMPERATURE  =  86  DEG  F  ?  _ 

RANGE  =  600  FT  ?83_ 

RANGE=83FT?_ 

IMPACT  VELOCITY  =  2800  FT/SEC  ?  3600_ 

IMPACT  VELOCITY  =  3600  FT/SEC  ?  _ 

SIGNAL  LEVEL  =  30  ?  60_ 

SIGNAL  LEVEL  =  60  ?_ 


"053194d3, 4500spm,  lOOrds,  prod  muz  rest,  h7-hl2’' 

OEHLER  SYSTEM  82  94-05-3 1  1 1 :45H  MACHINE  GUN  CHRONOGRAPH 


DISTANCE  =  10  MIN  VEL  =  3300  MAX  ROF  =  5000  TGT  TEMP  =  86 
FRAME  =  72.6  IN  SIGNAL  =  60  TGT  VEL  =  3600  RANGE  =  83 

1 1  ;45H  COMMENT:  053 194d3_ 

FIRE  BURST...THEN  TYPE  A  SPACE 

ROUND  EPOCH  A.TIME  VEL  RND/M  HORIZ  VERT  TIME2 


1 

0 

0 

0 

0 

-5.272 

-1,176 

0 

2 

0.02722 

0 

0 

2204 

-3.29 

0.716 

0 

3 

0.05019 

0 

0 

2611 

-5.512 

-1.21 

0 

4 

0.07153 

0 

0 

2811 

-3.175 

0.602 

0 

5 

0.09116 

0 

0 

3055 

-4.004 

1.817 

0 

6 

0.10893 

0 

0 

3377 

-5.267 

3.413 

0 

7 

0.12652 

0 

0 

3411 

-5.05 

1.422 

0 

8 

0.1428 

0 

0 

3686 

-5.244 

2.506 

0 

9 

0.15882 

0 

0 

3744 

-5.992 

-0.465 

0 

10 

0.17451 

0 

0 

3823 

-3.933 

0.839 

0 

11 

0.18966 

0 

0 

3959 

-3.911 

1.326 

0 

12 

0.20427 

0 

0 

4107 

-6.067 

2.08 

0 

13 

0.21899 

0 

0 

4076 

-5.897 

0.989 

0 

14 

0.23366 

0 

0 

4090 

-5.165 

1.807 

0 
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15 

0,24752 

0 

0 

4329 

-5.516 

-1.29 

0 

16 

0.26206 

0 

0 

4126 

-4.235 

-1.131 

0 

17 

0.27604 

0 

0 

4290 

-4.487 

0.864 

0 

18 

0.28977 

0 

0 

4371 

-5.506 

2.205 

0 

19 

0.30349 

0 

0 

4371 

-4.478 

1.194 

0 

20 

0.31702 

0 

0 

4433 

-3.97 

2.165 

0 

21 

0.33057 

0 

0 

4428 

-4.436 

1.177 

0 

22 

0.34402 

0 

0 

4460 

-3.089 

2.06 

0 

23 

0.35732 

0 

0 

4512 

-3.755 

3.244 

0 

24 

0.37052 

0 

0 

4544 

-4.418 

3.437 

0 

25 

0.38394 

0 

0 

4471 

-4.973 

1.727 

0 

26 

0.39723 

0 

0 

4514 

-3.752 

2,424 

0 

27 

0.41025 

0 

0 

4608 

-3.942 

0.576 

0 

28 

0.42332 

0 

0 

4592 

-2.741 

1.747 

0 

29 

0.43658 

0 

0 

4524 

-5.015 

1.727 

0 

30 

0.44972 

0 

0 

4565 

-5.903 

2.623 

0 

31 

0.46244 

0 

0 

4719 

-6.064 

1.543 

0 

32 

0.47552 

0 

0 

4587 

-5.641 

2.831 

0 

33 

0.48884 

0 

0 

4504 

-4.316 

1.262 

0 

34 

0.50171 

0 

0 

4660 

-3.711 

1.322 

0 

35 

0.51469 

0 

0 

4623 

-4.5 

1.3 

0 

36 

0.5273 

0 

0 

4756 

-5.981 

3.415 

0 

37 

0.54036 

0 

0 

4595 

-5.749 

2.42 

0 

38 

0.55341 

0 

0 

4594 

-5.692 

2.907 

0 

39 

0.56648 

0 

0 

4593 

-5.282 

0,059 

0 

40 

0.57944 

0 

0 

4626 

-3.767 

1.25 

0 

41 

0.59242 

0 

0 

4623 

-4.282 

1.542 

0 

42 

0.60521 

0 

0 

4692 

-6.261 

3.172 

0 

43 

0.61812 

0 

0 

4647 

-5.499 

1.544 

0 

44 

0.63114 

0 

0 

4608 

-5.268 

1.536 

0 

45 

0.64378 

0 

0 

4745 

-4.743 

0.406 

0 

46 

0.65687 

0 

0 

4584 

-4.622 

2.007 

0 

47 

0.66955 

0 

0 

4730 

-4.88 

0.635 

0 

48 

0.6829 

0 

0 

4495 

-4.895 

2.363 

0 

49 

0.6956 

0 

0 

4724 

-6.051 

1.48 

0 

50 

0.70858 

0 

0 

4621 

-4.942 

2.045 

0 

51 

0.72152 

0 

0 

4636 

-4.606 

1.761 

0 

52 

0.73416 

0 

0 

4747 

-4.451 

2.152 

0 

53 

0.74735 

0 

0 

4549 

-3,933 

0.635 

0 

54 

0.76055 

0 

0 

4546 

-5.469 

1.625 

0 

55 

0.77319 

0 

0 

4746 

-5.114 

1.155 

0 

56 

0.7862 

0 

0 

4609 

-5.067 

3.444 

0 

57 

0.79917 

0 

0 

4628 

-5.923 

-0.939 

0 

58 

0.81213 

0 

0 

4629 

-4.065 

1.41 

0 

59 

0.82514 

0 

0 

4610 

-4.457 

0.758 

0 

60 

0.83762 

0 

0 

4809 

-6.111 

3.05 

0 

61 

0.85074 

0 

0 

4570 

-5.738 

1.544 

0 

62 

0,8638 

0 

0 

4597 

-3.873 

1 

0 

63 

0.87656 

0 

0 

4702 

-4.756 

0.919 

0 

67 


64 

0.88965 

0 

0 

65 

0.90228 

0 

0 

66 

0.91539 

0 

0 

67 

0.92812 

0 

0 

68 

0,94097 

0 

0 

69 

0.95385 

0 

0 

70 

0.96685 

0 

0 

71 

0.97993 

0 

0 

72 

0.99282 

0 

0 

73 

1.00535 

0 

0 

74 

1.01849 

0 

0 

75 

1.03124 

0 

0 

76 

1.04429 

0 

0 

77 

1.05712 

0 

0 

78 

1.07 

0 

0 

79 

1.08283 

0 

0 

80 

1.09583 

0 

0 

81 

1.10876 

0 

0 

82 

1.12171 

0 

0 

83 

1.13452 

0 

0 

84 

1.14755 

0 

0 

85 

1.16036 

0 

0 

86 

1.17342 

0 

0 

87 

1.18634 

0 

0 

88 

1.19924 

0 

0 

89 

1.21225 

0 

0 

90 

1.22496 

0 

0 

91 

1.23794 

0 

0 

92 

1.2509 

0 

0 

93 

1.26402 

0 

0 

94 

1.27682 

0 

0 

95 

1.28986 

0 

0 

96 

1.30268 

0 

0 

97 

1.31576 

0 

0 

98 

1.32881 

0 

0 

99 

1.34154 

0 

0 

100 

1.3544 

0 

0 

101 

1.36708 

0 

0 

102 

1.37985 

0 

0 

103 

1.39286 

0 

0 

104 

1.40601 

0 

0 

105 

1.41905 

0 

0 

106 

1.43151 

0 

0 

107 

1.44458 

0 

0 

108 

1.45823 

0 

0 

109 

1.47066 

0 

0 

no 

1.48355 

0 

0 

111 

1.49666 

0 

0 

112 

1.50948 

0 

0 

4582 

-3.798 

2.946 

0 

4748 

-4.333 

1.05 

0 

4576 

-5.647 

2.259 

0 

4714 

-5.011 

2.002 

0 

4669 

-5.101 

1.693 

0 

4659 

-5.479 

2.527 

0 

4613 

-4.23 

1.838 

0 

4588 

-4.584 

1.647 

0 

4655 

-5.494 

3.357 

0 

4788 

-7.181 

1.195 

0 

4565 

-5.163 

2.748 

0 

4705 

-5,56 

0.317 

0 

4598 

-4.904 

1.651 

0 

4675 

-5.126 

1.905 

0 

4659 

-6.769 

0.672 

0 

4675 

-5.309 

2.086 

0 

4613 

-5.808 

0.867 

0 

4639 

-5.68 

0.3 

0 

4635 

6.474 

12.329 

0 

4684 

-4.201 

0.809 

0 

4604 

-5.763 

2.153 

0 

4682 

-5.344 

1.875 

0 

4594 

-5.562 

1.718 

0 

4641 

-5.67 

1.062 

0 

4654 

-3.758 

0.843 

0 

4611 

-4.824 

0.698 

0 

4719 

-5.78 

2.111 

0 

4624 

-4.065 

1.309 

0 

4627 

-4.673 

2.007 

0 

4575 

-4.837 

1.236 

0 

4684 

-3.821 

1.695 

0 

4601 

-4.509 

0.106 

0 

4681 

-5.096 

2.061 

0 

4585 

-6.078 

0.008 

0 

4600 

-5.341 

0.618 

0 

4711 

-4.671 

0.047 

0 

4665 

-4.358 

1.626 

0 

4734 

-4.841 

1.122 

0 

4698 

-5.696 

2.971 

0 

4611 

-5.52 

1.718 

0 

4563 

-6.452 

3.43 

0 

4598 

-5.094 

0.068 

0 

4815 

-7.604 

2.491 

0 

4592 

-6.462 

1.69 

0 

4395 

-5.264 

2.968 

0 

4827 

-5.88 

1.285 

0 

4653 

-6.023 

2.335 

0 

4577 

-5.572 

0.947 

0 

4681 

-5.187 

1.134 

0 
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113 

1.52248 

0 

0 

114 

1.53527 

0 

0 

115 

1.54833 

0 

0 

116 

1.56149 

0 

0 

117 

1.57391 

0 

0 

118 

1.58703 

0 

0 

119 

1.59959 

0 

0 

120 

1.61268 

0 

0 

121 

1.62558 

0 

0 

122 

1.63876 

0 

0 

123 

1.65169 

0 

0 

124 

1.66465 

0 

0 

125 

1.67718 

0 

0 

126 

1.69019 

0 

0 

127 

1.70343 

0 

0 

128 

1.71609 

0 

0 

129 

1.7287 

0 

0 

130 

1.7419 

0 

0 

131 

1.75461 

0 

0 

132 

1.76804 

0 

0 

133 

1.78088 

0 

0 

134 

1.79375 

0 

0 

135 

1.80654 

0 

0 

136 

1.81993 

0 

0 

137 

1.83273 

0 

0 

138 

1.84554 

0 

0 

139 

1.8585 

0 

0 

140 

1.87149 

0 

0 

141 

1.88452 

0 

0 

142 

1.89762 

0 

0 

143 

1.91051 

0 

0 

144 

1.9234 

0 

0 

145 

1.93645 

0 

0 

146 

1.94943 

0 

0 

147 

1.96305 

0 

0 

148 

1.97636 

0 

0 

149 

1.9901 

0 

0 

150 

2.00493 

0 

0 

151 

2.02053 

0 

0 

152 

2.03734 

0 

0 

153 

2.05527 

0 

0 

154 

2.07477 

0 

0 

155 

2.09557 

0 

0 

156 

2.11907 

0 

0 

157 

2.14536 

0 

0 

158 

2.17518 

0 

0 

159 

2.2044 

0 

0 

160 

4.3306 

0 

0 

161 

4.49671 

0 

0 

4614 

-5.263 

1.549 

0 

4689 

5.118 

13.656 

0 

4596 

-6.146 

0.558 

0 

4558 

-5.226 

0.927 

0 

4831 

-6.524 

-0.465 

0 

4571 

-4.168 

-0.555 

0 

4778 

-5.338 

2.582 

0 

4585 

-6.601 

1.031 

0 

4650 

-6.631 

0.765 

0 

4552 

-5.647 

2.107 

0 

4638 

-5.941 

0.101 

0 

4631 

-4.648 

1.905 

0 

4785 

-6.172 

0.642 

0 

4613 

-6.178 

2.106 

0 

4530 

-5.675 

0.059 

0 

4740 

-4.82 

2.986 

0 

4758 

-5.794 

1.624 

0 

4544 

-2.903 

-0.441 

0 

4722 

-4.911 

2.43 

0 

4465 

-6.352 

1.073 

0 

4672 

-6.018 

0.693 

0 

4664 

-4.706 

2.77 

0 

4691 

-4.596 

2.143 

0 

4480 

-4.234 

2.309 

0 

4687 

-4.638 

2.567 

0 

4684 

-5.72 

2.094 

0 

4629 

-4.487 

-0.178 

0 

4618 

-6.088 

2.33 

0 

4604 

-4.577 

0.377 

0 

4579 

-5.077 

0.068 

0 

4656 

-3.933 

0.957 

0 

4654 

-6.257 

1.496 

0 

4597 

-6.134 

0.156 

0 

4622 

-5.241 

2.036 

0 

4405 

-4.987 

-0.157 

0 

4505 

-4.359 

0.406 

0 

4369 

-4.858 

1.57 

0 

4045 

-5.432 

2.743 

0 

3845 

-6.76 

0.477 

0 

3567 

-5.466 

0.495 

0 

3347 

-3.137 

-0,462 

0 

3077 

-3.895 

0.978 

0 

2884 

-4.457 

-1.164 

0 

2553 

-5.683 

0.93 

0 

2282 

-3.801 

1.148 

0 

2011 

-3.129 

-0.182 

0 

2053 

-5.213 

-0,085 

0 

28 

20000 

20000 

0 

361 

20000 

20000 

0 
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SUMMARY:  RNDS  1  - 161  2134  RND/M 
MEAN  STD  DEV  MAX  MIN  RANGE 
VELOCITY  0  0  0  0  0 

HORIZONTAL  -4.925  1.512  6.474-7,604  14.078 
VERTICAL  1.52  1.684  13.656-1.29  14.946 
RADIAL  2.263 

RNDS  1  -  1594300 
RNDS  160  -  161  361 

"ENTER  SUB-BURST  START  AND  END  ROUND  NUMBERS  SEPARATED  BY  A  COMMA  12,8 1_" 


SUMMARY:  RNDS  12-  81  4577  RND/M 
MEAN  STD  DEV  MAX  MIN  RANGE 
VELOCITY  0  0  00  0 

HORIZONTAL  -4.975  0.854  -2.741  -7.181  4.44 
VERTICAL  1.632  1.015  3.444  -1.29  4.734 
RADIAL  1.327 

"ENTER  SUB-BURST  START  AND  END  ROUND  NUMBERS  SEPARATED  BY  A  COMMA  83,1 13_" 

SUMMARY:  RNDS  83-  113  4639  RND/M 
MEAN  STD  DEV  MAX  MIN  RANGE 
VELOCITY  0  0  0  0  0 

HORIZONTAL  -5.265  0.831  -3.758  -7.604  3.845 
VERTICAL  1.474  0.875  3.43  0.008  3.422 
RADIAL  1.207 

"ENTER  SUB-BURST  START  AND  END  ROUND  NUMBERS  SEPARATED  BY  A  COMMA  115,150_ 

SUMMARY:  RNDS  115  -  150  4599  RND/M 
MEAN  STD  DEV  MAX  MIN  RANGE 
VELOCITY  0  0  0  0  0 

HORIZONTAL  -5.306  0.874  -2.903  -6.631  3.728 
VERTICAL  1.217  1.077  2.986-0.555  3.541 
RADIAL  1.387 

ENTER  SUB-BURST  START  AND  END  ROUND  NUMBERS  SEPARATED  BY  A  COMMA  _ 

"H:45H  COMMENT:  dispersion  0.93,  2  flyersj' 
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APPENDIX  B.  CROSS  CORRELATION  CODE 


This  appendix  contains  the  cross  correleation  code  used  to  better  estimate  the  time 
of  the  first  round.  The  intent  was  to  use  the  XCORR.m  function  in  the  Signal  Processing 
Toolbox.  XCORR.m  cross  correlates  two  files  using  the  following  equation. 


N-\m\  „.o 

Due  to  the  size  of  the  files  there  were  memory  problems,  and  the  command 
function  XCORR.m  had  to  be  modified.  CROSSCORR.m,  is  a  modified  version  of  the 
function.  The  purpose  of  CROSSCORR.m  is  to  obtain  a  better  estimate  of  the  time  of  the 
first  round  fired.  CROSSCORR.m  uses  the  rough  estimated  from  THRESHOLD.m  as  a 
bench  mark  to  cross  correlate  the  times  between  the  accelerometers  and  dispersion.  For 
cross  correlation  two  files  are  used.  One  is  the  sum  of  the  squares  of  the  horizontal  and 
vertical  muzzle  restraint  accelerations.  The  second  contains  the  time  of  each  round  and  its 
horizontal  and  vertical  dispersion.  CROSSCORR.m  cross  correlates  for  three  hundred 
values  around  the  rough  estimate,  fifty  prior  and  two  hundred  and  fifty  after.  It  uses  the 
following  equation  to  cross  correlate 

-  40960- |m  I 

xcor= -  xcor(m,V)*x(n,l)*batacc(n*m,l)  (2) 

40,960- Iml  „.i 

where  m  is  the  value  for  the  bin  around  the  rough  estimate,  and  batacc  is  the  sum  of  the 
squares.  The  correlated  time  is  used  as  the  best  estimate  of  when  the  PHALANX  gun 
began  to  fire. 
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%  crosscorr.m 

%  Calculate  when  the  first  bullet  is  fired  by  using  the  xcorr  function  in  the  Signal  Proccesing  Toolbox 

clear  all; 
pit; 

load  batx 

xsense=9.73/1000;  %  V/g 

xscale=l/(xsense/386);  %  in/sec'^2  per  V 
batx=xscale*batx; 
tinc=l/10240; 

ttotal=(length(batx)- 1  )*tinc; 
tx=0:tinc:ttotal; 

load  baty; 

ysense=9.63/1000;  %  V/g 

yscale=l/(ysense/386);  %  in/sec^2  per  V 

baty=yscale*baty; 

tinc=l/10240; 

ttotal=(length(baty )- 1  )*tinc; 

ty=0:tinc:ttotal; 

t=tx; 

batacc=(batx.'^2  +  baty.''2);  %  Sum  of  the  Squares 
plot(t,batacc) 

title('Sum  of  the  Squares  (Muzzle  Restraint)’) 
xlabel(['Time,  sec  ’,num2str(date)]) 
ylabel(’Acceleration  ^^2,  (in/sec^2)'^2)’) 
grid 
pause 

load  firedata3; 
for  i=l:159 

firedata3(i,l)=firedata3(i,l)*10240+10117; 

end 

shott=zeros(40960,3); 

fori=l:159 

j=round(firedata3(i,  1)); 
shott(j,l)=firedata3(i,l); 
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shott(j  ,2)=firedata3  (i  ,2); 
shott(j,3)=firedata3(i,3); 

end 

x=zeros(size(batx)) ; 

for  i=l:40960 

if  abs(shott(i,2))  >  0; 

x(i)=l; 

end 

end 

plot(t,x) 

title(’Normalized  Dispersion  Displacement') 

xlabel([Time,  sec  ',num2str(date)]) 

ylabel('Normalized  Displacemnt*) 

grid 

pause 


clear  batx  batxl  batx2  batx3  batx4  batx5  batx6  batx7  batxS  batx9  batxlO 
clear  baty  baty  1  baty2  baty3  baty4  baty5  baty6  baty7  batyS  baty9  batylO 
clear  xsense  ysense  xscale  yscale  tx  ty  ttotal  firedata3  shott 

%  Estimate  the  cross-correlation  between  the  firadata3  and  the  sum  of  the  squares  of  the  muzzle  restraint  acceleration 
(vertical  and  horizontal) 

tinput=input('Enter  the  estimated  time  of  the  first  bullet  firing  '); 

bin=round(tinput/tinc); 

bin 

for  m=bin-50:bin+250 
m 

xcor(m,l)=0.0; 
for  n=l:40960-abs(m) 

xcor(m,  1  )=xcor(m,  1  )+x(n,  1  )*batacc(n+m,  1 ) ; 
end 

xcor(m,  1  )=xcor(m,  1  )/(40960-abs(m)); 
end 

xcor(bin-50:bin-i-250)*tinc; 

[maxcor,mmaxcor]=max(xcor) 
bullet=(mmaxcor+ 1 )/ 1 0240 
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plot(xcor(bin-50:bin+250)*tinc) 
titleCCross  Correlation  Between  FiredataS  &  batacc') 
xlabelCBin  Number') 
ylabelCCross  Correlation’) 

text(l  15,700,  ['First  Bullet  Time  ',num2str(bullet)]) 
grid 
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APPENDIX  C.  GTESTFFT.M 


This  appendix  contains  the  modified  code  written  by  Mike  Hatch.  GTESTFFT.M 
converts  the  acceleration  from  volts  to  meters/second''2.  It  also  computes  the  normalized 
FFT.  It  plots  the  acceleration  versus  time  and  the  normalized  FFT  versus  frequency. 


75 


%  gtestfftm  calculate  the  fft  of  forcing  function 
clear  all; 
pit; 

load  afty.mat; 

ysense=  10.25/1000;  %  V/g 

yscale  =  1/(0.102);  %  Tn/sec^2/g 

ysense=:yscale/ysense; 

afty=ysense*afty; 

tine  =  1/10240; 

ttotal  =  (iength(afty)-l)*tinc; 

ty  =  0:tinc:ttotal; 

plot(ty,afty) 

titleCAFT  BEARING  ACCELERATION  (VERTICAL)  ') 
xlabel([Time,  sec  -  Fire  3  ’,num2str(date)]) 
ylabel(Acceleration,  in/sec^2') 
grid 

f  =  10240*(1:20480)/40960; 
pn  =  abs(fft(afty))*2/length(afty); 
plot(f,pn( 1 :20480)) 

title(PSD  OF  AFT  BEARING  ACCELERATION  (VERTICAL)') 
xlabel(['Frequency,  Hz  ',num2str(date)]) 
ylabelC  Acceleration,  in/sec'^2’) 
grid 

%  load  and  plot  x  direction 
load  aftx.mat; 

xsense  =  10.24/1000;  %  W/g 

xscale=  1/(0.102);  %  m/sec'^2/g 

xsense=xscale/xsense; 

aftx  =  xsense*aftx; 

tine  =  1/10240; 

ttotal  =  (length(aftx)-l)*tinc; 

tx  =  0:tinc:ttotal; 

plot(tx,aftx) 

tideCAFT  BEARING  ACCELERATION  (HORIZONTAL)  ’) 
xlabel([Time,  sec  -  Fire  3  ’,num2str(date)]) 
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ylabelC  Acceleration,  in/sec^2') 
grid 

f  =  10240*(l:20480)/40960; 
pn  =  abs(fft(aftx))*2/length(aftx); 
plot(f,pn( 1:20480)) 

titleCPSD  OF  AFT  BEARING  ACCELERATION  (HORIZONTAL)’) 
xlabel(['Frequency,  Hz  ',num2str(date)]) 
ylabel('Acceleration,  in/sec^2’) 
grid 

%  load  and  plot  X  direction 
load  batx.mat; 

xsense  =  9.73/1000;  %  V/g 

xscale=  1/(0.102);  %  mfsec^lfg 

xsense=:xscale/xsense; 

batx  =  xsense*batx; 

tine  =  1/10240; 

ttotal  =  (length(batx)-l)*tinc; 

tx  =  0:tinc:ttotal; 

plot(tx,batx) 

titleCMUZZLE  RESTRAINT  ACCELERATION  (HORIZONTAL)  ’) 
xlabel([Time,  sec  -  Fire  3  ',num2str(date)]) 
ylabel('AcceIeration,  in/sec^2') 
grid 

f  =  10240*(l:20480)/40960; 
pn  =  abs(fft(batx))*2/length(batx); 
plot(f,pn(l:20480)) 

dtleCPSD  OF  MUZZLE  RESTRAINT  ACCELERATION  (HORIZONTAL)’) 
xlabel(['Frequency,  Hz  ',num2str(date)]) 
ylabel('Acceleration,  in/sec'^2’) 
grid 


%  load  and  plot  y  direction 
load  baty.mat; 

ysense  =  9,63/1000;  %  V/g 
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yscale=  1/(0.102);  %  m/sec^2/g 

ysense=yscale/ysense; 

baty  =  ysense*baty; 

tine  =  1/10240; 

ttotal  =  (length(baty)-l)*tinc; 

ty  =  0:tinc:ttotal; 

plot(ty,baty) 

titleCMUZZLE  RESTRAINT  ACCELERATION  (VERTICAL) ') 
xlabel([Time,  sec  -  Fire  3  ',num2str(date)]) 
ylabel(Acceleration,  in/sec^2') 
grid 

f  =  10240*(l:20480)/40960; 
pn  =  abs(fft(baty))*2/length(baty); 
plot(f,pn(l:20480)) 

titleCPSD  OF  MUZZLE  RESTRAINT  ACCELERATION  (VERTICAL)') 
xlabel([Frequency,  Hz  ',num2str(date)]) 
ylabel(Acceleration,  in/sec'^2’) 
grid 

%  load  and  plot  X  direction 
load  midx.mat; 

xsense  =  10,09/1000;  %  V/g 

xscale=  1/(0.102);  %  ni/sec^2/g 

xsense=xscale/xsense; 

midx  =  xsense*midx; 

tine  =  1/10240; 

ttotal  =  (iength(midx)- 1  )*tinc; 

tx  =  0:tinc:ttotal; 

plot(tx,midx) 

titleCMID  HOUSING  ACCELERATION  (HORIZONTAL)  ’) 
xlabel([Time,  sec  -  Fire  3  ',num2str(date)]) 
ylabel(Acceleration,  in/sec'^2') 
grid 

f  =  10240*(1:20480)/40960; 

pn  =  abs(fft(midx))*2/length(midx); 
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plot(f,pn(  1:20480)) 

titleCPSD  OF  MID  HOUSING  ACCELERATION  (HORIZONTAL)') 
xlabel(['Frequency,  Hz  ’,num2str(date)]) 
ylabel('Acceleration,  in/sec^2') 
grid 

%  load  and  plot  y  direction 
load  midy.mat; 

ysense  =  10.02/1000;  %  V/g 

yscale=  1/(0.102);  %  m/sec^2/g 

ysense=yscale/ysense; 

midy  =  ysense*midy; 

tine  =  1/10240; 

ttotal  =  (length(niidy)”l)*tinc; 

ty  =  0:tinc:ttotai; 

plot(ty,midy) 

titleCMID  HOUSING  ACCELERATION  (VERTICAL)  ’) 
xlabel([Time,  sec  -  Fire  3  ',nuin2str(date)]) 
ylabelC  Acceleration,  in/sec^2') 
grid 

f  =  10240*(1:20480)/40960; 

pn  =  abs(fft(midy))*2/length(inidy); 

plot(f,pn(l:20480)) 

titleCPSD  OF  MID  HOUSING  ACCELERATION  (VERTICAL)’) 
xlabel(['Frequency,  Hz  ',num2str(date)]) 
ylabelC  Acceleration,  in/sec'^2’) 
grid 

%  load  and  plot  z  direction 
load  midz.mat; 

zsense=  10.05/1000;  %  V/g 

zscale=  1/(0.102);  %  mlscc^Hg 

zsense=zscale/zsense; 

midz  =  zsense*niidz; 

tine  =  1/10240; 

ttotal  =  (length(inidz)-l)*tinc; 

tz  =  0:tinc:ttotal; 
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plot(t2,midz) 

titieCMID  HOUSING  ACCELERATION  (AXIAL) ') 
xlabel([Time,  sec  -  Fire  3  ',num2str(date)]) 
ylabel(Acceleration,  in/sec^2') 
grid 


f  =  10240*(l:20480)/40960; 

pn  =  abs(fft(midz))*2/length(midz); 

plot(f,pn(l:20480)) 

title(PSD  OF  MID  HOUSING  ACCELERATION  (AXIAL)’) 
xlabel([’Frequency,  Hz  ’,num2str(date)]) 
ylabel('Acceleration,  in/sec^2') 
grid 

%  load  and  plot  signal 
load  sig.mat; 
tine  =  1/10240; 
ttotal  =  (length(sig)-l)*tinc; 
ts  =  0:tinc:ttotal; 

plot(ts,sig) 

titleCnRING  SIGNAL  ’) 

xlabel([Time,  sec  -  Fire  3  ’,num2str(date)]) 

ylabelCAcceleration,  in/sec'^2') 

grid 
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APPENDIX  D.  THRESHOLD.M 


This  appendix  contains  the  MATLAB  code,  THRESHOLD. m.  It  sets  a  limit 
value  to  determined  the  first  occurrence  of  an  acceleration  magnitude.  THRESHOLD.m 
returns  the  bin  number  when  the  limit  is  exceeded.  This  gives  a  rough  estimate  of  the 
time  of  the  first  round  fired.  This  estimate  is  used  to  calculate  a  better  initial  firing  time 
with  CROSSCORR.m. 
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%  Determine  when  the  first  bullet  is  fired. 
%  threshold.m 
clear ; 
pit; 

load  batx3.txt; 

threshold=0.25; 

for  i=l:  12288 

if  batx3(i)  >  abs(threshold) 

bx=i; 

break 

end 

end 

load  baty3.txt; 

fori=l:12288 

if  baty3(i)>  abs(threshold) 

by=i; 

break 

end 

end 

load  aftx3.txt; 

for  i=l:  12288 

if  aftx3(i)>  abs(threshold) 

ax=i; 

break 

end 

end 

load  afty3.txt; 

for  i=l:  12288 

if  afty3(i)>  abs(threshold) 

ay=i; 

break 

end 

end 

load  midx3.txt; 

for  i=l:  12288 

if  mi<ix3(i)>  abs(threshold) 

mx=i; 

break 

end 
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end 

load  midy3.txt; 

for  i=l:  12288 

if  midy3(i)>  abs(threshold) 

my=i; 

break 

end 

end 

load  inidz3.txt; 

for  i=l:  12288 

if  niidz3(i)>  abs(threshold) 

mz=i; 

break 

end 

end 
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APPENDIX  E.  CALIBRATION  CERTIFICATES  FOR  THE  QUARTZ  SHEAR 
MODE  ICP  ACCELEROMETERS  USED  FOR  DATA  COLLECTION  ON  THE 
PHALANX  GUN  SYSTEM. 


This  appendix  contains  the  calibration  charts  for  the  Series  353  accelerometers. 
Table  E-1  summarizes  the  accelerometer  set  up  for  31  May  1994. 


Accelerometer 

Model 

Serial  Number 

Location 

Direction 

r353B03 

SN  5966 

Mid  Housing 

Axial 

1353B03 

SN  5967 

Mid  Housing 

Vertical 

1353B03 

SN  5968 

Mid  Housing 

Horizontal 

1353B03 

SN9319 

Aft  Bearing 

Horizontal 

1353B03 

SN  5969 

Aft  Bearing 

Vertical 

353B03 

SN11851 

Muzzle  Restraint 

Horizontal 

SN  11850 

Vertical 

Table  E-1,  31  May  1994  Accelerometer  Set  Up 


85 


Piezotronics,  Inc.  3425  Walden  Avenue  Depew,  NY  14043-2495  USA 


oc=i  « 

tr 

opr!  <f 
<=5  CO 

^  to 

0 


Q  <M 

Wo'? 

^  O  ^ 

o  II 

9  -v  II 
CC  w  .. 

tr.  E  ^ 


0^0)^ 

+! 


CO 

in 

z 

cu 

o 

+ 

s 

H 

CO 

G5 

in 

< 

s 

• 

CO 

o 

in 

CO 

1 

LL 

CD 

o 

UJ 

Q. 

CO 

o 

c 

o 

B 

CD 

C 

03 

QC 

>* 

O) 

c 

o 

0) 

cL 

E 

UJ 

05 

CD 

CD 

CC 

CC 

h- 

N 

X  . 

^  (J>  > 


s  S) 

G3  ro 
in 


s  ^ 

Q  I  • 

S) 

m 


S  CO 

Q 

m 


Q  « 
ro  ,• 


cu 

G)  CO  CO 

•  •  •  • 

O  Q  ^  o  C3D 


< 

>> 

K 

>s 

< 

Q 

;> 

U) 

o 

c: 

CD 

CD 

> 

c 

13 

> 

Z 

o 

U) 

CD 

CO 

cr 

CD 

V.- 

c 

03 

CD 

-J 

H 

CD 

CD 

U- 

to 

CO 

03 

in 

B 

< 

GC 

m 

CO 

CD 

CD 

CO 

CD 

> 

CO 

c 

CO 

c 

c 

o 

o 

< 

B 

O 

c 

2 

o 

CO 

CD 

CD 

E 

CL 

o 

> 

CC 

H 

o 

Piezotronics,  Inc.  3425  Walden  Avenue  Depew.  NY  14043-2495  USA 


Piezotronics,  Inc.  3425  Waiden  Avenue  Depew,  NY  14043-2495  USA 

716-684-0001  Calibrated  by 


Piezotronics,  Inc.  3425  Waiden  Avenue  Depew,  NY  14043-2495  USA 

716-684-0001 


Piezotronics,  Inc.  3425  Walden  Avenue  Depew,  NY  14043-2495  USA 


Piezotronics,  Inc.  3425  Walden  Avenue  Depew,  NY  14043-2495  USA 

716-684-0001 


APPENDIX  F.  TIME  DOMAIN  WINDOW 


This  appendix  contains  the  MATLAB  code  used  to  apply  a  time  domain  window 
to  the  acceleration.  The  function  uses  a  normalized  half  cosine  and  half  sine  and 
multiplies  it  by  the  acceleration  in  the  time  domain. 
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function  y  =  timwin(tl,t2,t3,t4,timdat,xdat); 

%  timwin.m 

%  Create  a  window  that  when  applied  to  a  data  it  will  set  zeros  up 
%  to  an  initial  time,  then  apply  half  a  cosine  and  increase  the 
%  value  from  zero  to  one  then  stay  at  one  until  a  certain  time 
%  then  decrease  back  to  zero  with  half  a  cosine  wave. 

y=zeros(size(xdat)); 
ndx  1  =timndx(t  1  ,timdat); 
ndx2=timndx(t2,timdat) ; 
ndx3=timndx(t3,timdat); 
ndx4=timndx(t4,timdat) ; 

y(ndx2:ndx3,l)  =  ones(ndx3-ndx2+l,l); 

denom=ndx2-ndx  1 ; 
for  i=ndxl:ndx2 

y(i,  1  )=0.5*(1  -cos(pi*  (i-ndx  1  )/denom)) ; 
end 

denom=ndx4-ndx3 ; 
for  i=ndx3:ndx4 

y(i,  1  )=0.5*(  1  +cos(pi*(i’ndx3)/denom)); 
end 


y  =  y  .*  xdat; 
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APPENDIX  G.  FREQUENCY  DOMAIN  WINDOW 


This  appendix  contains  the  MATLAB  code  used  to  apply  a  frequency  domain 
window  to  the  FFT.  The  function  uses  a  normalized  half  cosine  and  half  sine  and 
multiplies  it  by  the  FFT  in  the  frequency  domain. 
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function  y  =  freqwin(fl,f2,fdat,fftdat); 

%  freqwin.m 

%  Create  a  window  that  when  applied  to  fft  data  it  will  set  zeros 
%  up  to  an  initial  time,  then  apply  half  a  cosine  and  increase  the 
%  value  from  zero  to  one  then  stay  at  one  until  a  certain  time 
%  then  decrease  back  to  zero  with  half  a  cosine  wave. 

y=zeros(size(fftdat)); 

lenfdat=length(fdat); 

finc=0.25; 

ndxl=round(f  l/fmc)+l ; 
ndx2=round(f2/finc)+l ; 
ndx3=lenfdat-ndx2+2; 
ndx4=lenfdat-ndx  1  +2 ; 


y(ndx2:ndx3,l)  =  ones(ndx3-ndx2+l,l); 

denom=ndx2-ndx  1 ; 
for  i=ndxl:ndx2 

y(i,  1  )=0.5*(  1  “COs(pi*(i‘ndx  1  )/denom)); 
end 

denom=ndx4-ndx3 ; 
for  i=ndx3:ndx4 

y(i,  1  )=0.5*(  1  +cos(pi*(i-ndx3)/denom)); 
end 

y  =  y  .*  fftdat; 
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APPENDIX  H.  CUBIC  SMOOTHING  SPLINES  CODE 


CSAPSXY.m  is  code  that  applies  the  CSAPS.m  command  of  the  Splines 
Toolbox,  to  the  muzzle  restraint  displacement  curves  for  the  horizontal  and  vertical 
directions.  It  plots  the  displacement  curve  obtained  from  direct  integration  along  with  the 
cubic  spline  approximation. 
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%  Calculate  the  Cubic  Smoothing  of  the  Muzzle  Restraint  Displacement 
%  (Horizontal)  obtained  from  the  integration  algorithim. 

%  Using  MatLab  csaps  function 
%  Repeat  for  a  third  p  value 

clear  all; 
pit; 

load  time; 
load  disintx; 

tO=cputime; 

%  Plot  the  Muzzle  Restraint  Displacement  (Horizontal)  vs  Time 

figure(l) 

plot(time,disx) 

title(’Muzzle  Restraint  Displacement  (Horizontal)’) 
xlabel([Time,  sec  -  Fire  3  ’,num2str(date)]) 
ylabel('Displacement,  m') 
grid 
hold  on 

%  Calculate  and  plot  the  Cubic  Spline  for  the  third  p  value 
%  Plot  and  compare  to  the  Displacement  curve  obtained  from  the  integration 
%  Increase  the  Cubic  Spline  Parameter  by  1000000000 

epsilon=max(diff(time))^3/16; 
p3x=l/(l+epsilon*  1000000000); 
p3x 

fdisxl2=csaps(time,disx,p3x,time); 
ddisx  1 2=disx-fdisx  1 2; 
save  fdisxl2  fdisxl2; 
save  ddisx  12  ddisx  12; 

%  Plot  the  Cubic  Smoothing  Spline  and  print 
plot(time,fdisx  1 2,’r') 

text(  1,1.1,  'Cubic  Spline  Parameter  epsilon*  10000000') 
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plot(time,dciisxl  2,'g') 
print 

figure  (2) 
plot(tinie,ddisxl2) 

title('Muzzle  Restraint  Displacement  (Horizontal)') 
xlabel([Time,  sec  -  Fire  3  ',num2str(date)]) 
ylabel('Displacement,  m') 
grid 

el-(cputime-t0)/3600 

%  Calculate  the  Cubic  Smoothing  of  the  Muzzle  Restraint  Displacement 
%  (Vertical)  obtained  from  the  integration  algorithim. 

%  Using  MatLab  csaps  function 
%  Repeat  for  a  third  p  value 

t3=cputime; 

load  disinty; 

%  Plot  the  Muzzle  Restraint  Displacement  (Vertical)  vs  Timeload  disinty; 

figure  (3) 

plot(time,disy) 

title('Muzzle  Restraint  Displacement  (Vertical)’) 
xlabel([’Time,  sec  -  Fire  3  ‘,num2str(date)]) 
ylabelCDisplacement,  m’) 
grid 
hold  on 

%  Calculate  and  plot  the  Cubic  Spline  for  the  third  p  value 
%  Plot  and  compare  to  the  Displacement  curve  obtained  from  the  integration 
%  Increase  the  Cubic  Spline  Parameter  by  1000000000 

epsilon=max(diff(time))'^3/ 1 6 ; 
p3y=l/(l+epsilon*  1 000000000); 
p3y 

fdisyl  2=csaps(time,disy,p3y,time): 
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ddisy  1 2=disy-fdisy  1 2; 
savefdisyl2fdisyl2; 
save  ddisy  12  ddisy  12; 

%  Plot  the  Cubic  Smoothing  Spline  and  print 
plot(time,fdisy  1 2,  'r') 

text(  1,2.1,  'Cubic  Spline  Parameter  epsilon*  10000000’) 
plot(time, ddisy  1 2,’g*) 


figure  (4) 
plot(time, ddisy  12) 

title('Muzzle  Restraint  Displacement  (Vertical)’) 
xlabeI([’Time,  sec  -  Fire  3  ',num2str(date)]) 
ylabel('Displacement,  m') 
grid 

e3=(cputime-t3)/3 600 
ttotai=(cputime-tO) 
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